ASD-TDR-63-405 


STUDY  TO  FIND  SIMULANTS  FOR  FUELS 
FOR  USE  IN  STRUCTURES  FATIGUE  TESTING 


TECHNICAL  DOCUMENTARY  REPORT  No.  ASD-TDR-63-405 


NOVEMBER  1963 


AIR  FORCE  FLIGHT  DYNAMICS  LABORATORY 
RESEARCH  AND  TECHNOLOGY  DIVISION 
AIR  FORCE  SYSTEMS  COMMAND 
WRIGHT-PATTERSON  AIR  FORCE  BASE,  OHIO 


Project  1347,  Task  134704 


(Prepared  under  Contract  AF  33(657)-7901  by 
Rocketdyne,  Canoga  Park,  California.) 

^oo(,omo03 


NOTICES 


When  Government  drawings,  specifications,  or  other  data  are  used  for  any 
purpose  other  than  in  connection  with  a  definitely  related  Government  procure¬ 
ment  operation,  the  United  States  Government  thereby  incurs  no  responsibility 
nor  any  obligation  whatsoever;  and  the  fact  that  the  Government  may  have 
formulated,  furnished,  or  in  any  way  supplied  the  said  drawings,  specifications, 
or  other  data,  is  not  to  be  regarded  by  implication  or  otherwise  as  in  any 
manner  licensing  the  holder  or  any  other  person  or  corporation,  or  conveying 
any  rights  or  permission  to  manufacture,  use,  or  sell  any  patented  invention 
that  may  in  any  way  be  related  thereto. 


Qualified  requesters  may  obtain  copies  of  this  report  from  the  Defense 
Documentation  Center  (DDC),  (formerly  ASTIA),  Cameron  Station,  Bldg.  5, 
5010  Duke  Street,  Alexandria  4,  Virginia 


This  report  has  been  released  to  the  Office  of  Technical  Services,  U.S. 
Department  of  Commerce,  Washington  25,  D.C.,  in  stock  quantities  for  sale 
to  the  general  public. 


Copies  of  this  report  should  not  be  returned  to  the  Aeronautical  Systems 
Division  unless  return  is  required  by  security  considerations,  contractual 
obligations,  or  notice  on  a  specific  document. 


200  -  January  1964  -  162-22-396 


FOREWORD 


This  report  was  prepared  by  the  Rocketdyne  Division  of 
North  American  Aviation,  Inc.,  under  USAF  Contract  No. 
AF33(657)-7901.  This  contract  was  initiated  under 
Project  Nr.  1347,  "Structural  Testing  of  Flight  Vehicles," 
Task  134704,  "Thermal  Application  and  Control."  The  work 
was  administered  under  the  direction  of  the  Air  Force 
Flight  Dynamics  Laboratory,  Research  and  Technology 
Division,  with  Mr.  Clark  E.  Beck  acting  as  Project 
Engineer. 


This  report  covers  work  conducted  from  April  1962  to 
January  1963. 


Contributors  to  this  report  were  W.  Bisbee,  J.  E.  Flana¬ 
gan,  L.  Grant,  C.  Hamermesh,  J.  Hilzinger,  C.  Fujikawa, 
and  E.  A.  Lawton. 


ABSTRACT 
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INTRODUCTION 


When  a  vehicle  is  in  flight,  its  structure  is  externally  heated  by  con¬ 
vection  from  the  boundary  layer.  However,  due  in  part  to  various  internal 
heat  sinks,  the  structural  temperatures  will  not  be  uniform.  Therefore, 
the  temperature  gradients  and  resultant  differential  thermal  expansions 
will  produce  thermal  stresses  in  the  structure.  Such  environments  or 
stresses  must  be  reproduced  during  a  test  program  if  the  structural  integ¬ 
rity  of  the  vehicle  is  to  be  established  conclusively.  The  thermal  effects 
of  the  fuel  or  propellant  mass  must  be  considered  as  a  heat  sink,  and 
accounted  for  during  structural  testing.  However,  the  use  of  actual  fuel 
or  propellant  during  a  full-scale  structural  test  would  be  dangerous,  both 
to  the  operators  conducting  the  tests  and  the  equipment  employed.  This 
equipment  is  extensive,  expensive,  and  not  readily  moved  to  a  location 
where  live  propellants  might  be  employed.  Therefore,  simulants  for  live 
propellants,  both  of  the  liquid  and  solid  type,  are  required  for  vehicle- 
structure  verification  testing.  The  major  criteria  for  satisfactory  pro¬ 
pellant  simulation  are:  (l)  safety  in  use,  as  reflected  by  little  or  no 
fire  or  toxicity  hazard,  and  (2)  duplication  of  propellant  properties. 


Manuscript  released  by  the  author  22  April  1963  for  publication  as  an  ASD 
Technical  Documentary  Report. 
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SUMMARY 


A  program  was  performed  to  obtain  nonhazardous  simulants  which  dupli¬ 
cated  the  pertinent  physical  properties  of  several  liquid  and  solid  pro¬ 
pellants  of  current  interest.  Materials  which  appear  to  simulate  best  the 
five  liquid  propellants  studied  are: 

1.  For  monomethylhydrazine ;  water,  dibutylcarbitol ,  diethylcarbitol , 
and  butylcellosolve 

2.  For  the  Titan  II  fuel  blend;  water,  dibutylcarbitol,  diethylcarbitol, 
and  butylcellosolve 

3.  For  JP-4;  dibutylcarbitol 

4.  For  RP-1 ;  dibutylcarbitol  and  red  oil  (commercial  oleic  acid) 

5.  For  nitrogen  tetroxide;  Freon  -21  and  Freon  -11 

Other  materials  which  might  serve  as  simulants  under  special  conditions  are 
also  listed. 

Simulation  of  an  aluminized  solid  propellant  containing  ammonium  perchlor¬ 
ate  and  a  polybutadiene  binder  was  achieved  by  replacement  of  the  perchlor¬ 
ate  in  the  coarse  oxidizer  phase  by  ammonium  sulfate.  The  simulant  was 
completely  nonhazardous  and  duplicated  the  mechanical  and  thermal  proper¬ 
ties  of  the  "live”  propellant  satisfactorily. 

Thermal  properties,  particularly  the  diffusivity  of  solid  propellants,  are 
generally  not  measured.  Therefore,  an  evaluation  of  the  test  methods  used 
for  nonpropellant  application  was  made,  and  nonhazardous  techniques  useful 
for  propellant  application  were  devised. 
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Calculation  of  the  specific  heat  of  a  propellant  from  the  specific  heats  of 
its  ingredients  is  feasible,  and  represents  a  nonhazardous  means  for  obtain¬ 
ing  this  parameter. 

Attempts  to  use  salts  such  as  potassium  chloride  or  sodium  chloride  in  lieu 
of  ammonium  sulfate  in  the  simulant  propellant  resulted  either  in  no  cure 
or  in  only  a  partial  cure.  This  was  attributed  to  the  greater  absorption 
of  the  crosslinker  by  the  two  metal  salts. 
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SIMULATION  OF  LIQUID  PROPELLANTS 


By  mutual  agreement  between  representatives  of  the  Air  Force  and  Rocket- 
dyne,  the  following  liquids  were  selected  for  simulation:  HP-1,  mono- 
methylhydrazine ,  the  Titan  II  fuel  blend,  nitrogen  tetroxide,  and  the  jet 
fuel  JP-4 . 

Complete  liquid-phase  physical  properties  data  for  these  propellants  were 
not  available.  Therefore,  the  liquid  simulant  phase  of  this  program  was 
divided  into  three  parts:  (i)  a  literature  search  to  compile  all  of  the 
known  data  on  the  selected  propellants;  (2)  an  evaluation  of  techniques 
for  correlating  and  extrapolating  pertinent  physical  properties;  and  (3) 
a  selection  of  simulants  and,  where  required,  an  estimation  of  their 
properties. 

For  the  propellant  selected,  those  properties  and  the  temperature  range 
over  which  they  must  be  duplicated  by  the  simulant  will  depend  to  a  large 
extent  on  the  type  of  heating  the  propellant  is  likely  to  encounter.  This 
heating  may  be  of  two  types;  (l)  equilibrium,  and  (2)  transient.  Equi¬ 
librium  heating  refers  to  heating  which  is  generally  slow;  for  example, 
the  heating  or  cooling  of  propellants  in  the  normal  environmental  tem¬ 
perature  range  of  -65  to  +200  F.  Heating  or  cooling  does  not  occur 
rapidly  but  is  accomplished  over  a  relatively  long  period  of  time.  In 
direct  contrast  is  transient  heating.  Here,  the  skin  of  the  vehicle  is 
heated  rapidly;  i.e.,  large  thermal  gradients  occur.  For  the  latter  type 
of  heating,  temperatures  as  high  as  500  F  might  be  considered  for  effec¬ 
tive  testing  of  the  vehicle  structure. 
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The  properties  that  appear  most  important  from  the  standpoint  of  equilibrium 
heating  are  specific  heat,  viscosity,  coefficient  of  thermal  expansion,  and 
vapor  pressure.  For  the  transient  type  of  heating,  properties  such  as  ther¬ 
mal  conductivity,  heat  of  evaporation,  and  critical  state  become  important 
and  must  be  simulated. 

To  provide  the  best  simulant  for  each  of  the  selected  propellants  over  the 
applicable  temperature  ranges,  the  properties  of  the  propellants  in  these 
ranges  must  be  known.  A  comprehensive  search  of  all  available  literature 
for  the  known  applicable  physical  properties  of  the  selected  propellants 
was  therefore  made.  The  results  of  this  survey  are  reported  in  the  perti¬ 
nent  sections,  and  presented  as  the  change  with  temperature  of  the  specific 
propellant  property.  In  most  cases,  such  properties  have  been  determined 
or  estimated  for  near-ambient  conditions  only.  Thus,  the  data  available 
fell  short  of  the  temperature  ranges  of  interest  to  this  study.  This  situ¬ 
ation  had  been  anticipated;  therefore,  concurrent  with  the  literature  survey, 
an  evaluation  of  various  empirical  or  semiempirical  techniques  for  estimat¬ 
ing  physical  properties  was  made.  Extrapolation  of  available  data  was  made 
based  on  what  was  believed  to  be  the  best  technique  available. 


CORRELATION  OF  PHYSICAL  PROPERTIES 

The  prediction  of  physical  properties  is  usually  based  on  the  correlation 
of  known  information,  with  interpolation  or  extrapolation  employed  as 
required.  Correlations  may  be  divided  into  three  types:  wholly  empirical, 
semiempirical  (i.e.,  based  on  some  theoretical  concept),  and  wholly  theo¬ 
retical.  The  latter  type,  although  most  desirable,  is  usually  lacking  for 
liquids  because  of  the  difficulties  inherent  in  the  understanding  of  the 
liquid  state.  The  purely  empirical  correlation  suffers  from  unreliability, 
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particularly  when  it  has  been  developed  for  one  class  of  compounds  only* 
Most  of  the  useful  correlations  fall  into  the  semiempirical  class;  i.e., 
they  are  of  a  form  suggested  in  part  by  theory,  with  empirical  constants 
based  on  experimental  data. 

The  value  of  a  method  of  estimation  depends  on  the  accuracy,  simplicity, 
and  the  type  of  information  required  for  its  use.  If  sufficient  informa¬ 
tion  is  available,  it  may  be  possible  to  estimate  the  property  quite 
accurately.  On  the  other  hand,  if  there  is  little  information,  a  less 
reliable  estimate  is  obtained. 

A  measure  of  the  degree  of  confidence  in  the  accuracy  of  such  estimation 
techniques  for  the  properties  studied  was  obtained  from  a  comparison  of 
the  experimental  and  calculated  values  in  the  temperature  range  where 
such  experimental  data  existed.  No  detailed  statistical  analysis  of  the 
comparison  was  made,  but  simple  average  errors  (without  regard  to  sign) 
are  indicated  where  applicable. 


Critical  State  Constants 


Many  estimation  techniques  are  based  on  the  correlation  of  various  phys¬ 
ical  properties  with  the  critical  pressure  and  temperature.  The  avail¬ 
ability  of  methods  for  either  calculating  or  determining  experimentally 
the  critical  pressures  and  temperatures  is  therefore  of  importance.  The 
methods  that  have  been  used  for  predicting  critical  temperatures  and  pres¬ 
sures  are,  to  a  large  degree,  empirical,  the  relationships  being  devised  to 
fit  the  experimental  data.  Such  an  approach  causes  some  uncertainty  in  the 
prediction  of  the  critical  constants  of  materials  which  are  different  in 
chemical  composition  from  those  for  which  the  techniques  were  developed. 
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In  those  cases  where  the  critical  constants  of  pure  materials  were  cal¬ 
culated,  the  method  of  Lydersen  (Ref.  1  )  was  used.  The  method  of  Lyder- 
sen  and  that  of  Riedel  (Ref.  2  )  are  of  comparable  accuracy  (Ref.  3), 
both  methods  being  based  on  the  sunmation  of  atomic  and  structural  con¬ 
stants  which  represent  the  component  parts  of  the  molecule. 

For  mixtures  such  as  the  Titan  II  fuel  blend  and  for  some  of  the  simulants, 
Kay’s  rule  (Ref.  4  )  appears  to  be  most  appropriate  for  calculation  of  the 
pseudocritical  pressures  and  temperatures.  Kay’s  rule  is  a  simple,  or 
linear  mixture  rule  which  assumes  that  the  critical  constant  of  the  mix¬ 
ture  is  equal  to  the  sum  of  the  products  of  mole  fractions  and  the  criti¬ 
cal  constants  of  the  pure  components. 


Vapor  Pressure 


The  pressure  exerted  by  a  vapor  in  equilibrium  with  the  liquid  is  known  as 
the  vapor  pressure.  Because  the  equilibrium  between  a  liquid  and  its  vapor 
is  dependent  on  temperature,  equations  relating  vapor  pressure  and  tempera¬ 
ture  are  desirable.  In  this  respect,  the  correlation  of  Riedel  (Ref.  5  ) 
is  useful  because  it  applies  to  a  variety  of  materials  over  a  wide  range 
of  temperature.  The  Riedel  equation  is 


log10  it  “  0  (Tr)  +  <«c  -  7)  *  (Tr)  (1) 

where  J i  (T^)  and  lb  (T^)  are  functions  of  the  reduced  temperature.  Deri¬ 
vation  of  Riedel’s  equation  for  a  specific  material  requires  experimental 
determination  of  the  vapor  pressure  at  a  single  temperature.  From  this 
one  value,  Of  may  be  determined  from  Eq.  1,  since  Riedel  has  tabulated 
(T^)  and  (T^)  as  functions  of  T  . 
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Reidel's  equation  has  been  applied  to  various  substances  over  the  entire 
liquid  range  from  the  freezing  to  the  critical  point.  When  tested  by  use 
of  the  smoothed  experimental  data  reported  by  Stull  (Bef.  6),  the  calcu¬ 
lated  and  experimentally  determined  vapor  pressures  agreed  vithin  ±10  per¬ 
cent  for  most  compounds.  This  agreement  is  good,  considering  that  the 
pressures  ranged  from  1  ran  Hg  to  values  greater  than  100,000  mm  Hg. 


Density 


At  least  one  value  of  the  orthobaric  density  is  usually  available  for 
almost  any  compound,  thus,  Watson's  density  expansion  factor  (Bef.  7  and 
8)  can  be  used.  Watson  found  that  the  quotient  of  density  divided  by 
the  expansion  factor,  <4  was  essentially  constant  for  all  liquids.  The 
expansion  factor,  when  expressed  as  a  function  of  the  reduced  temperature 
and  pressure,  is  approximately  the  same  for  all  substances.  Thus,  evalu¬ 
ation  of  the  expansion  factor  as  a  function  of  the  reduced  pressure  and 
temperature  for  one  substance  for  which  there  is  sufficient  data,  permits 
the  use  of  Eq.  2  for  predicting  the  densities  of  other  compounds. 

- 

P  -  w  u  <2> 

1 

All  that  is  required  to  predict  the  density  of  any  compound  by  this 
equation  is  one  density  value  at  one  tenqierature ,  so  that  can  be 

evaluated.  Using  this  method,  the  values  usually  are  well  within  ±5  per¬ 
cent,  even  for  polar  substances  such  as  water.  The  greatest  deviations 
usually  occur  in  the  critical  region. 
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Viscosity 


A  large  number  of  viscosity-temperature  relations  have  been  proposed  and 
tested.  However,  the  method  of  Uyehara  and  Watson  (Ref,  9)  appears  to 
be  the  most  useful  because  it  is  applicable  to  many  types  of  compounds 
as  well  as  mixtures.  This  correlation  is  based  on  the  reduced  viscosity 
(t?  -  7)/v  )  and  the  reduced  temperature  and  pressure.  By  use  of  the  gen¬ 

eralized  plot  of  Uyehara  and  Watson,  the  critical  viscosity  of  a  material 
may  be  obtained  if  a  single  experimental  viscosity  is  known.  For  those 
cases  where  no  viscosity  data  are  available,  useful  approximations  of  vis¬ 
cosity  may  be  made. 

While  no  entirely  satisfactory  method  is  available  for  the  estimation  of 
viscosities  of  mixtures,  Uyehara  and  Watson's  method  shows  fair  agreement 
with  experimental  data  for  gaseous  mixtures.  It  is  believed  that  this 
method  is  simpler  and  more  reliable  than  those  methods  employing  the  vis¬ 
cosities  of  pure  components  at  the  conditions  of  the  mixture. 

The  basic  equation  by  which  viscosities  may  be  estimated  as  a  function  of 
temperature  is 

*  =  VrVc  (3) 


Thermal  Conductivity 

A  knowledge  of  the  thermal  conductivity  is  important  in  transport  phe¬ 
nomena.  Attempts  to  derive  thermal  conductivity  from  purely  theoretical 
grounds  have  achieved  limited  success.  Few  reliable  results  on  experi¬ 
mental  thermal  conductivity  measurements  were  available  until  recently. 
Furthermore,  most  of  the  reported  measurements  cover  only  a  narrow  tem¬ 
perature  range. 
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The  methods  available  that  involve  the  use  of  the  most  readily  accessible 
data  are  the  equations  of  Weber  (Ref.  10  )  and  Smith  (Ref.  11  ).  They  dif¬ 
fer  only  in  the  constant  of  Eq.  4,  and  are  basically  empirical. 


K 


£  1/3 

M 


where 

K  = 

the  thermal  conductivity 

C 

P 

the  specific  heat 

/>  - 

the  density 

M  = 

the  molecular  weight 

A  = 

a  constant 

(4) 


A  third  equation  is  available.  This  equation  differs  from  those  previ¬ 
ously  mentioned  in  that  Cp  is  the  specific  heat  at  Tc/2,  and  A  is  a  "uni¬ 
versal  constant"  for  unassociated  liquids  (Lapushkin,  Ref. 12  ).  If  the 
product  ACpM  (Eq.  5  )  is  experimentally  determined, 

K  =  ACpM-l/3(^)4/3  (5) 


it  can  be  set  equal  to  a  constant  B,  and  results  in 


K  = 


(6) 


A  comparison  of  calculated  results  (Eq.  6  )  with  experimental  data  for 
liquid  hydrocarbons  and  associated  liquids  (water  and  alcohols)  gave 
deviations  less  than  5  percent  (Ref. 12  ).  Larger  deviations  are  obtained 
for  mixtures. 
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While  the  method  of  Sakiadis  and  Coates  (Ref.  13)  has  a  firmer  theo¬ 
retical  basis,  the  reliability  of  the  results  obtained  are  not  much  bet¬ 
ter  than  those  from  the  more  empirical  equations.  Furthermore,  this 
method  involves  the  use  of  data  which  are  not  always  available  nor 
easily  estimated. 

The  Cragoe  equation  appears  best  for  hydrocarbons  (Ref.  14),  and  was  used 
for  RP-1  and  JP-4.  For  all  other  systems,  the  equation  employed  was  that 
which  gave  the  closest  agreement  for  the  temperature  range  of  existing 
experimental  data.  For  those  substances  for  which  no  experimental  data 
existed,  the  La pushkin  equation  was  used  in  most  cases.  The  reliability 
of  the  results  however,  is  no  better  than  the  accuracy  of  the  density  and 
specific  heat  values  used. 


Heat  Capacity 

Owing  to  the  difficulties  in  dealing  with  the  liquid  state  in  general, 
little  progress  has  been  made  toward  the  prediction  of  heat  capacities 
of  liquids  on  a  theoretical  basis.  Those  empirical  correlations  which 
are  reported  in  the  literature  are  often  restricted  to  specific  classes 
of  compounds  and  sometimes  require  knowledge  of  physical  property  data 
which  may  not  be  available. 

The  heat  capacities  of  liquid  petroleum  fractions  have  been  studied,  and 
correlations  reported  (Ref.  14).  However,  because  most  of  the  compounds 
and  mixtures  which  were  investigated  under  this  contract  differed  in  con¬ 
stitution  from  petroleum  fractions,  such  a  specific  method  could  not  be 
used  with  confidence.  A  more  general  method  which  might  be  more  appli¬ 
cable  to  the  series  of  compounds  under  consideration  was  sought. 
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The  method  of  Chow  and  Bright  (Ref.  15)  correlates  the  variation  of 
with  temperature  by  use  of  Watson1  s  density  expansion  coefficient*  The 
correlation  suggested  by  Chow  and  Bright  is 

C  CO  a  =  b  (7) 

P 

where  the  exponent  "a"  was  assigned  a  value  of  2.8  on  the  basis  of  a 

study  of  some  representative  alcohols,  aromatic  hydrocarbons,  and  halides. 

The  term  "b"  varies  with  the  nature  of  the  substance,  but  within  a  given 

homologous  series,  is  reasonably  constant.  The  use  of  Chow  and  Bright’s 

correlation  for  the  present  application  involved  evaluating  "a"  and  "b" 

for  each  substance  studied.  This  was  possible  by  determining  the  slope 

"a"  of  the  straight  line  obtained  from  a  plot  of  log  7r““  vs  log  tOf  the 

P 

intercept  nb,!  was  then  calculated  from  the  equation  of  the  straight  line. 
The  determination  of  "a"  and  "bM  permitted  the  calculation  of  the  depend¬ 
ence  of  C  on  T.  Good  agreement  between  experimental  and  calculated  data 
was  observed  in  the  temperature  range  for  which  experimental  data  existed. 


Miscellaneous  Physical  Properties 


Properties  such  as  freezing  points  and  boiling  points  are  usually  avail¬ 
able.  Heats  of  vaporization  at  the  boiling  point  can  be  calculated  from 
a  knowledge  of  the  normal  boiling  point  (Ref.  7  and  16  ).  Knowledge  of 
the  variation  of  density  with  temperature  permits  calculation  of  the  coef¬ 
ficient  of  thermal  expansion  by  the  relationship; 


where 

p, 

p. 


the  t^-t^  density  at  t^ 
the  density  at  t 
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PROPELLANT  PHYSICOCHEMICAL  PROPERTIES 


The  following  sections  summarize  the  physicochemical  properties  (as 
experimentally  determined  or  from  the  correlations  discussed  above)  of 
monomethylhydrazine,  the  Titan  II  fuel  blend,  RP-1,  JP-4,  and  nitrogen 
tetroxide . 


MONOMETHYLHYDRAZINE 

Commercial  monomethylhydrazine  (MMH)  is  95+  percent  CH  NHNH  .  The  usual 

j  £ 

impurities  are  methylamine  and  water.  MMH  is  a  clear,  water-white, 
hygroscopic  liquid  which  tends  to  turn  yellow  upon  exposure  to  air.  It 
is  a  toxic,  volatile  liquid  which  absorbs  CO^  from  the  air  at  ambient 
temperatures.  It  is  completely  miscible  in  all  proportions  with  hydra¬ 
zine,  low  molecular  weight  alcohols,  and  hydrocarbons.  It  is  not  sensi¬ 
tive  to  impact  or  friction.  MMH  is  more  stable  than  hydrazine  on  mild 
heating,  and  similar  to  hydrazine  in  its  sensitivity  to  catalytic 
oxidation. 


Physicochemical  Properties 

A  summary  of  the  physicochemical  properties  of  MMH  is  given  in  Table  1. 


Critical  Constants.  The  critical  constants  of  MMH  have  been  experiment¬ 
ally  determined  by  the  Aerojet-General  Corporation  (Ref.  17)-  The  criti¬ 
cal  constants  are  T  =  594  F,  P  -  1195  psia,  and  P  =  18.1  lb/  cu  ft. 

c  c  c 
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TABLE  1 


SUMMARY  OF  PHYSICO-CHEMICAL  PROPERTIES  OF  MONOMETHYLHYDRAZINE 


Property 

Value 

Reference 

Molecular  formula 

CH^NIINHg 

Molecular  weight 

46.075 

Melting  point,  F 

-62.5 

18 

Boiling  point,  F 

189 

18 

Flash  point,  F 

54 

19 

Density,  liquid,  at  77  F,  Ib/cu  ft 

54.  58 

17,  20 

Viscosity,  liquid,  at  77  F,  lb/ft-sec 

51.8  x  10"5 

17,  20 

Vapor  pressure  at  77  F,  psia 

0.96 

18,  21 

Critical  temperature,  F 

594 

17 

Critical  pressure,  psia 

1195 

17 

Heat  of  vaporization  at  77  F,  Btu/lb 

376.9 

18 

Heat  of  formation,  liquid,  at  77  F,  Btu/lb 

496 

22 

Specific  heat,  liquid,  at  77  F,  Btu/lb-F 

0.700 

18 

Thermal  conductivity,  Btu/hr-ft-F, (Calc) 

0.16 

12 

16 


Vapor  Pressure.  The  vapor  pressures  of  MMH  over  the  temperature  range  32 
to  77  F  have  been  determined  by  Aston  (Ref.  18  );  Aerojet  (Ref.  17)  extended 
these  measurements  to  the  critical  point.  The  results  (Fig.  1  )  are  satis¬ 
factorily  represented  (average  deviation,  2.9  percent  from  1  to  1195  psia) 
by  Eq.  9  . 

log  jr  =  0  (Tr)  +  0.290  0  (Tr)  (9) 

r 

Density.  The  density  of  MMH  in  the  temperature  range  -62  to  +181  F, 
under  its  own  vapor  pressure,  has  been  experimentally  determined  (Ref.  17). 
Extrapolation  of  these  data  to  the  critical  point  was  made  by  use  of  the 
equation: 

^  (lb/cu  ft)  =  414.2  CO  (l0) 

The  average  deviations  of  the  calculated  values  from  the  experimentally 
determined  densities  were  within  1  percent.  The  calculated  critical  den¬ 
sity  of  18.2  lb/cu  ft  agrees  quite  favorably  with  the  experimental  value 
of  18.1  lb/cu  ft.  Because  of  this  agreement  it  is  felt  that  the  equation 
above  will  give  good  estimates  of  the  density  of  MMH  throughout  the  liq¬ 
uid  range,  under  its  own  vapor  pressure  (Fig.  2). 

Viscosity.  Both  the  Aerojet-General  and  Metalectro  Corporation  (Ref.  17 

and  20)  measured  the  viscosity  of  MMH  as  a  function  of  temperature  over 

the  range  -62  to  +176  F.  Extension  of  these  data  to  the  critical  point 
was  made  by  use  of  Eq.  11. 

T)  (lb/ft-sec)  =  2.43  x  10  ^  TJr  (ll) 
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Figure  1.  Vapor  Pressure  of  Monomethy lhydraz ine 
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In  the  temperature  range  32  to  17 6  F,  the  average  deviation  between  cal¬ 
culated  and  experimental  data  is  6  percent.  The  absolute  and  kinematic 
viscosities  are  shown  in  Fig.  3  - 

The  accuracy  of  the  derived  equation  was  checked  at  210  F  and  250  F. 

The  kinematic  viscosity  was  measured  at  210  F  and  found  to  be  0.378  cen- 
tistokes;  at  250  F,  the  measured  value  was  0.310  centistokes.  The  devi¬ 
ation  from  the  extrapolated  value  at  210  F  is  3  percent,  and  the  measured 
value  at  250  F  is  the  same  as  that  calculated. 


Thermal  Conductivity.  No  experimental  thermal  conductivity  measurements 
on  MMH  have  been  made.  The  thermal  conductivities  given  in  Fig.  4  were 
calculated  by  the  method  of  Lapushkin. 


Heat  Capacity.  The  heat  capacity  of  liquid  MMH  from  its  melting  point  to 
77  F  was  measured  by  Aston  (Ref.  18).  Extrapolation  of  the  experimental 
data  to  the  critical  point  (Fig.  5  )  was  made  by  use  of  Eq.  12. 

cp  aP'554  =  0.340  (12) 

The  average  deviation  between  experimental  and  calculated  values  in  the 
temperature  range  -45  to  +77  F  is  0.4  percent. 


Coefficient  of  Thermal  Expansion.  The  coefficients  of  thermal  expansion 
in  the  temperature  range  -50  F  to  the  critical  point  were  calculated  by 
means  of  Eq.  8.  The  results  are  shown  in  Table  2.  In  the  temperature 
range  60  to  79  F,  the  coefficient  of  thermal  expansion  has  been  reported 
as  -5- 918  x  10_4/F  (Ref.  23). 
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Figure  3-  Viscosity  of  Liquid  MMH 
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KINEMATIC  VISCOSITY  (CENT 1ST OKES) 


TABLE  2 


COEFFICIENTS  OF  THERMAL  EXPANSION  FOR  MONOMETHYLHYDRAZINE 


Temperature  Range,  F 

-P  x  104/F 

“62  to  -5.8 

5.29 

-5.8  to  77.0 

6.71 

77.0  to  180. 5 

6.58 

200  to  300 

8.74 

300  to  400 

11.1 

400  to  500 

11.5 

300  to  594 

106 
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TITAN  II  FUEL  BLEND 


The  50/50  blend  of  unsymmetrical  dimethylhydrazine  (UDMH)  and  hydrazine, 
a  mixture  of  toxic  materials,  has  been  selected  as  a  fuel  for  the  Titan 
II  ballistic  missile.  The  trade  name  of  Aerozine-50  has  been  adopted  by 
Aerojet-General  Corp.,  and  is  governed  by  the  Military  Specification 
MIL-P-27402  (USAF),  dated  25  August  1961.  The  propellant  specification 
is  given  in  Table  3  (Ref.  24). 


Physicochemical  Properties 

A  summary  of  the  physicochemical  properties  of  the  Titan  II  fuel  blend 
is  presented  in  Table  4 . 


Critical  State  Constants.  The  psuedocritical  state  constants  of  the 

50/50  fuel  blend  were  calculated  by  Kay*s  rule  (Ref.  4  ).  Values  used 

for  the  critical  constants  of  N_H.  and  UDMH  were  those  available  in  the 

2  4 

literature  (Ref.  21,  25,  and  26  ). 


Vapor  Pressure.  The  vapor  pressure  of  a  mixture  such  as  the  50/50  fuel 
blend  varies  with  the  ullage  of  the  system.  Experimental  data  are  avail¬ 
able  in  the  temperature  range  14  to  160  F  (Ref. 27  and  28  )  for  an  ullage 
of  46  percent.  Extrapolation  of  the  vapor  pressures  under  these  condi¬ 
tions  was  made  to  the  pseudocritical  point  by  means  of  Eq.  13. 

log  =  0  (Tr)  -  0.556  0  (Tr)  (13) 

r 

The  average  deviation  between  the  experimental  and  calculated  data  in  the 
temperature  range  32  to  160  F  was  3-5  percent.  The  published  and  extrapo¬ 
lated  data  are  shown  in  Fig.  6. 
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TABLE  3 

PROPELLANT  SPECIFICATION 
MIL-P-27402 


Chemical  Requirements 

Specification  Vt .  Percent 

N2H4 

51.0  ±0.8 

UDMH  plus  amines 

47 • 0  (min . ) 

H^O  plus  other  soluble 
impurities 

1.8  (max.) 

N0H  -  UDMH  plus 

98.2  (min.) 

^  1  amines 
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TABLE  4 


SUMMARY  OF  PHYSICO-CHEMICAL  PROPERTIES  OF  THE 
TITAN  II  FUEL  BLEND  (Ref.  24) 


Property 

Typical  Value 

Melting  point  range0,  F 

17.6  to  21.2 

Boiling  point  (incipient)^,  F 

149  F 

Density  at  14.7  psia  and  77  Fa,  lb/cu  ft 

56.1 

Viscosity  at  77  Fa,  lb/ft-sec 

54.9  x  10~5 

Vapor  pressure  at  77  Fc  (46  percent  ullage),  psia 

2.75 

Pseudocritical  temperature  (Calc),  F 

634 

Pseudocritical  pressure  (Calc),  psia 

1696 

Heat  of  vaporization  (Calc),  Btu/lb 

425.8 

Heat  of  formation  at  77  F  (Calc),  Btu/lb 

507.35 

Specific  heat  at  77  F,  Btu/lb-F 

0.684 

Thermal  conductivity  at  150  F^,  Btu/ft-hr-F 

0.165 

Measured  on  samples  of  the  fuel  blend  of  typical  composition 

(51.0  wt.  percent  NgH^,  wt.  percent  UDMH,  and  0.5  wt  percent  Ho0). 

^The  fuel  blend  is  not  a  constant  boiling  mixture. 

£ 

Fuel  blend  composition  51.0  wt.  percent  N0H. ,  48.4  wt.  percent  UDMH,  and 
0.6  wt.  percent  H^O.  ’ 

^Rocketdyne-determined  value. 
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Density.  The  density  of  the  50/50  fuel  blend  was  measured  by  the  Aerojet- 
General  Corp.  (Ref.  27)  to  160  F.  Extrapolation  of  these  data  to  550  F 
(Fig.  7  )  was  carried  out  by  means  of  Eq.  14  (Fig.  7). 

p  (lb/cu  ft)  =  420.2  to  (14) 


Viscosity.  The  viscosity  of  the  50/50  fuel  blend  was  measured  by  Aerojet- 
General  (Ref.  27  )  to  160  F.  Extrapolation  of  these  data  to  500  F  was 
made  by  means  of  Eq.  15 . 

T J  (lb/ft-sec)  =  2.15  x  10  5  ^  (15) 

Additional  measurements  were  made  by  Jtocketdyne  under  this  contract  at 
210  and  250  F.  For  the  50/50  fuel  blend,  the  measured  values  were  0.394 
centistokes  at  210  F,  and  0.330  centistokes  at  250  F.  These  data  dif¬ 
fered  from  the  extrapolated  values  by  -9  percent  at  210  F  and  -12  percent 
at  250  F  (Fig.  8). 


Thermal  Conductivity.  No  thermal  conductivity  measurements  of  the  50/50 
fuel  blend  have  been  reported.  During  the  course  of  this  contract,  ther¬ 
mal  conductivity  measurements  were  made  at  Rocketdyne.  The  following 
results  were  obtained: 


t(f) 

K(Btu/ft-hr-F) 

113.54 

0.163 

113.55 

0.167 

149.51 

0.165 

149.47 

0.164 
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,  lb/cu  ft 
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Figure  7.  Orthobaric  Density  of  Titan  II  Fuel  Blend 
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KINEMATIC  VISCOSITY  (CENT 1ST OKES) 


The  constant  B  of  the  Lapushkin  equation  was  determined  from  3*..  w  by 
substitution.  The  derived  equation  for  the  thermal  conductivity  of  the 
fuel  blend  is: 

K(Btu/ft-hr-F)  -  0.202  p (l6l 

The  extrapolated  data  are  shown  in  Fig.  9  • 


S :.,^ific  Heat.  Beil  has  measured  the  specific  heat  of  the  i\iz~  blent:  ar 
81.2  and  101.8  F  (lief.  28  )  by  the  method  of  mixtures.  From  the  measure., 
values  of  0.692  and  0.698  Btu/lb-F  and  the  Chow-Bright  correic oi , 

Eq.  17  was  derived  for  the  variation  of  heat  capacity  with  temperature 
The  data  are  plotted  in  F.g.  10  . 


C  CO0-690 
P 


0.172 


07; 


Coefficient  of  Thermal  Expansion.  The  coefficient  of  thermal  ex.mr.sion 
was  calculated  from  Eq.  8  .  The  results  are  shown  in  Table  5. 


RP-1 


BP-1  is  a  straight-run  kerosene  fraction  which  has  been  treated  to  remove 
aromatics  and  sulfur-containing  hydrocarbons.  It  consists  principally  of 
napthenic  and  paraffinic  hydrocarbons.  BP-1  is  generally  a  clear,  water- 
white  liquid.  It  is  chemically  stable  and  insensitive  to  shock.  It  wall 
react  with  strong  oxidizers  and,  at  conditions  of  extreme  pressure  and 
temperature ,  also  undergo  pyrolysis.  The  fuel  is  flammable  and  its  vapors 
form  explosive  mixtures  with  air.  Spills  of  rocket  oxidizers  and  BP-1  form 
mixtures  that  can  be  exploded  by  mechanical  shock,  heat,  or  spark. 
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1.10 


o 


Figure  10.  Specific  Heat  of  Liquid  Aerozine-50 
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TABLE  5 


COEFFICIENTS  OF  THERMAL  EXPANSION  OF  THE 
TITAN  II  FUEL  BLEND 


■  —  . - — . . — 1 

Temperature  Range,  F 

-fi  x  104/F 

20  to  100 

5.8 

100  to  200 

6.1 

200  to  300 

8.1 

| 

300  to  400 

8.7 

400  to  500 

n.9 

500  to  550 

18.0 
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Physicochemical  Properties 


A  summary  of  the  physicochemical  properties  of  RP-1  is  given  in  Table  6. 


Critical  State  Constants.  Since  RP-1  is  a  mixture  of  hydrocarbons,  the 
techniques  of  estimating  the  pseudocritical  properties  must  be  used.  In 
addition  to  Kay  (Ref.  4  ),  Smith  and  Watson  (Ref.  29),  and  Edmister 
(Ref.  30)  have  developed  methods  for  graphically  predicting  the  pseudo¬ 
critical  properties  of  hydrocarbon  mixtures.  The  critical  constants  of 
RP-1  (Table  6 )  are  the  average  values  obtained  by  use  of  these  three 
methods . 


Vapor  Pressure.  From  the  typical  vapor  pressure  of  0.3  psia  at  160  F 
for  RP-1,  Eq. 18was  derived. 

log  ~  =  <f>  (Tr)  =  0.043  tb  (Tr)  (18) 

r 

This  equation,  which  is  simpler  to  manipulate  than  that  derived  from  the 
Reid  vapor  pressure  equation  (Ref.  31  )  gives  comparable  calculated 
results  (Fig.  ll). 

Density.  Experimental  densities  (Ref.  3l)  in  the  temperature  range  40 
to  100  F  have  been  determined.  From  these  values,  Eq.  19  was  derived. 

P(lb/cu  ft)  =  361.2  60  (19) 

In  the  experimental  temperature  range,  the  average  deviation  of  the  cal¬ 
culated  values  was  0.1  percent  (Fig*  12). 
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TABLE  6 


SUMMARY  OF  PHYSICOCHEMICAL  PROPERTIES  OF  HP- 1 

(Ref.  32) 


Property 

Typical  Value 

Symbol 

Molecular  weight 

175  (n  =  12.5) 

Melting  point,  F 

-55 

Distillation 

10  percent,  F 

377 

End  point ,  F 

507 

Flash  point,  F,  average 

139 

Specific  gravity,  average,  at  68  F 

0.806 

Density,  Ib/cu  ft,  at  68  F 

50.3 

Viscosity  at  60  F,  lb/ft-sec 

1.4  x  10~5 

Vapor  pressure  at  160  F,  psia 

0.3 

Pseudocritical  temperature,  F 

780 

Pseudocritical  pressure,  psia 

310 

Heat  of  vaporization  at  N.B.P.,  Btu/lb 

106 

Heat  of  formation,  Btu/(CH  n  ). 

1 .  "5  'b 

-10.375 

Specific  heat  at  68  F,  Btu/lb-F 

0.45 

Thermal  conductivity  at  68  F,  Btu/hr-ft-F 

0.083 
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Vi  sc os i ty .  The  critical  viscosity  of  HP-1  was  calculated  from  the  fol¬ 
lowing  equation: 


V 


c 


.7-7  M1/2  2/3 

T  l/6  c 
pc 


where 

Tj  =  critical  viscosity  in  micropoise 
T  psuedocritical  temperature,  K 

P  ~  liquid  density,  g  per  cc 


CO  =  expansion  factor 

M  =  molecular  weight,  g  per  mol 


(20) 


The  value  of  r/  ,  calculated  for  EP-1  is  0.0261  centipoises,  or  0.75  x 
10  J  lb/ft  -sec.  In  addition  to  the  previous  reported  experimental  vis¬ 
cosities  of  RP-1  in  the  temperature  range  40  to  100  F  (Ref.  3l)  >  two 
additional  data  points  for  the  kinematic  viscosity  were  obtained  at  210 
and  250  F.  The  measured  values  of  0.784  centistokes  at  210  F  and  0.642 
centistokes  at  250  F  differed  by  1  and  4  percent,  respectively,  from 
those  calculated  by  Eq.  21. 

rj  (  lb/ft-sec)  =  1.75  x  10"5  «  (21) 

The  results  obtained  from  Eq.2l  are  given  in  Fig.  13. 


Thermal  Conductivity.  The  calculated  thermal  conductivity  of  liquid 
RP-1  at  14.7  psia  is  shown  in  Fig.  14.  The  equation  of  Cragoe  (Ref.  14) 
developed  for  petroleum  liquids,  was  used  for  the  calculation. 
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,  lb/ft-sec  X  10 


TIC  VISCOSITY,  CENTISTOKES 


JP-4 

JP-4  is  one  of  a  number  of  jet  engine  fuels.  Basically,  the  fuel  is  com¬ 
posed  of  hydrocarbons  with  a  restriction  on  the  allowable  percentage  of 
aromatics  present.  A  range  of  the  physical  properties  that  may  be  encoun¬ 
tered  in  various  samples  of  JP^-4  is  given  in  Table  7. 

Since  JP-4  was  chosen  to  be  simulated  from  -60to+230  F,  its  physical  prop¬ 
erties  were  determined  only  for  this  temperature  range. 
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TABLE  7 


RANGE  OF  PHYSICAL  PROPERTIES  OF  JP-4 
(Ref.  34) 


Property 


Typical  Value 


Average  Molecular  Weight 
Viscosity  at  100  F,  centistokes 
Pseudocritical  temperature,  F 
Pseudocritical  pressure,  psia 
Specific  heat  at  100  F,  Btu/lb-F 
Thermal  conductivity  at  100  F,  Btu/ft-hr-F 
Heat  of  combustion,  Btu/lb 


110-160 
0.67  to  1.40 
575  to  710 
310  to  510 

0.475  to  0.524 
0.0788 

18,500  to  18,725 


45 


Physicochemica 1  Properties 


A  summary  of  the  physicochemical  properties  of  JP-4  is  presented  in 
Table  8. 


Critical  Constants.  The  pseudocritical  temperature  and  pressure  are 
those  reported  in  Ref. 34  • 


Vapor  Pressure .  The  average  true  vapor  pressure  of  JP-4  below  zero  F 
is  insignificant.  The  average  vapor  pressures  as  reported  (Ref.  33) 
in  the  temperature  range  Oto  230  F  can  be  satisfactorily  calculated 
from  Eq.  24  (Fig.  16). 

log  jT"  =  t  (Tr)  -  1.804  (Tr)  (24) 

r 


Density.  The  average  density  of  JP-4  as  reported  (Ref.  33)  is  repre¬ 
sented  by  Eq.  25-  The  density  is  plotted  as  a  function  of  temperature 
in  Fig. 17  . 

p  (lb/cu  ft)  =  351.3  W  (25) 

Viscosity.  The  temperature  range  of  interest  precludes  the  use  of  the 
Uyehara-Watson  correlation  for  the  estimation  of  the  viscosity  of  JP-4. 
Therefore,  the  average  viscosity  values  given  in  Ref.  34  were  used  for 
estimation  purposes  (Fig.  18). 
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TABLE  8 


SUMMARY  OF  PHYSICOCHEMICAL  PROPERTIES  OF  JP-4 
(Ref,  33  and  34) 


Property 

Typical  Value 

Symbol 

Melting  point,  F 

Distillation 

(fWn 
<  -76 

10  percent,  F 

144 

End  Point,  F 

487 

Specific  gravity,  at  60  F 

0.773 

Viscosity,  at  60  F,  lb/ft-sec 

6  x  10~4 

Pseudocritical  temperature,  F 

720 

Pseudocritical  pressure,  psia 

370 

Heat  of  vaporization,  at  300  F,  Btu/lb 

126 

Heat  of  combustion,  Btu/lb 

18,680 

Specific  heat,  at  60  F,  Btu/lb-F 

0.472 

Thermal  conductivity,  at  60  F,  Btu/hr-ft-F 

0.082 
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PRESSURE,  PSIA 


TEMPERATURE,  F 


Figure  l6.  Vapor  Pressure  of  JP-4 
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20  40  60  80 
LATURE,  F 

•baric  Density  of  JP-4 
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Kinematic  Viscosity,  Centistokes 


Thermal  Conductivity.  The  thermal  conductivity  of  JP-4  is  shown  as  a 
function  of  temperature  in  Fig.  19.  The  data  were  calculated  by  means 
of  Eq . 22  . 


Specific  Heat.  The  specific  heat  of  JP-4  as  calculated  from  Eq.  23  is 
given  in  Fig. 20  for  the  temperature  range  -60  to  +80  F. 


NITROGEN  TETROXIDE 


Nitrogen  tetroxide  is  a  very  reactive  and  toxic  oxidizer.  It  is  insen¬ 
sitive  to  mechanical  shock,  heat,  or  detonation.  Although  it  is  non¬ 
flammable  with  air,  it  can  support  combustion  with  combustible  materials. 


Commercially  available  nitrogen  tetroxide  consists  primarily  of  the 
tetroxide  in  equilibrium  with  a  small  amount  of  nitrogen  dioxide 
(Eq.  27). 


N2°4<- 


2N0_ 


(27) 


In  the  solid  state,  N^O^  is  colorless,  diamagnetic,  and  exists  appar¬ 
ently  entirely  as  the  tetroxide.  In  the  liquid  and  gaseous  states,  a 
brown  color  and  paramagnetism  appear.  These  are  intensified  as  the 
temperature  increases.  In  addition,  gas  density  measurements  show  that 
the  percentage  dissociation  increases  from  20  at  17  F  to  90  at  212  F.  The 
percent  of  N0^  present  in  NgO^  is  also  a  function  of  pressure  (Ref.  35  )• 
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As  a  result  of  the  dissociation  phenomenon  and  the  various  factors  which 
affect  it,  the  survey  pertaining  to  the  physical  properties  of  N^O^  has 
been  limited  to  those  observed  under  its  own  vapor  pressure.  Even  under 
these  conditions,  the  correlation  techniques  could  not  be  satisfactorily 
applied  in  the  case  of  the  experimental  vapor  pressure,  viscosity,  and 
thermal  conductivity  data.  However,  since  the  physical  properties  in 
most  cases  are  known  from  the  freezing  point  to  temperatures  approaching 
the  critical  temperature,  the  failure  to  successfully  correlate  these 
properties  in  equation  form  did  not  cause  a  serious  problem. 


Physicochemical  Properties 

A  summary  of  the  pertinent  physical  properties  is  given  in  Table  9.* 


Vapor  Pressure.  The  vapor  pressure  of  NQ0^  has  been  experimentally 
determined  over  the  temperature  range  12  to  300  F  (Ref. 35,37*  and  4l). 
These  results  are  shown  in  Fig.  21. 

Density .  The  density  of  liquid  N^O^  has  been  measured  from  32  to  300  F 
under  its  own  vapor  pressure  (Ref.  38and  43)*  These  results  are  con¬ 
sistent  with  Eq.  28,  with  a  maximum  deviation  of  5  percent  and  an  average 
deviation  of  2  percent  (Fig.  22). 

P  (lb/ cu  ft)  =  775-5  <*)  (28) 
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TABLE  9 


SUMMARY  OF  PHYSIC 0CHEMICAL  PROPERTIES 

OF  N204 

Property 

Typical 

Value 

Ref erenc 

Molecular  formula 

N2°4 

Molecular  weight 

92.016 

Melting  point,  F 

11.8 

35 

Boiling  point,  F 

70.1 

35 

Density  at  14.7  psia,  and  68  F,  lb/cu  ft 

90.54 

38 

Viscosity  of  liquid  at  77  F,  lb/ft-sec 

2.80  x  10"4 

37 

Vapor  pressure  at  68  F,  psia 

13.92 

35 

Critical  temperature,  F 

316.8 

36 

Critical  pressure,  psia 

1469 

37 

Critical  density,  lb/cu  ft 

34.77 

38 

Heat  of  vaporization  at  B.P.,  Btu/lb 

178 

35 

Heat  of  formation  at  77  F  (liquid),  Btu/lb-Mol 

-12,240 

40 

Specific  beat  at  77  F,  Btu/lb-F 

0.374 

41 

Thermal  conductivity  at  77  F  (the  bubble  point), 
Btu/ft-hr-F 

0.0755 

42 
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Viscosity.  The  viscosity  of  N^O^  in  the  liquid  phase  has  been  measured 
as  a  function  of  temperature  from  14  to  280  F  (Ref. 39  and  43).  It  was 
not  possible,  however,  to  correlate  these  data  by  any  of  the  standard 
methods.  The  data  were  nevertheless  extrapolated  to  the  critical  tem¬ 
perature  by  means  of  Eq.  3,  using  the  experimental  value  at  280  F  as 
the  reference  viscosity  (Fig.  23). 

Thermal  Conductivity.  The  thermal  conductivity  of  liquid  N^O^  (Fig.  24) 
in  the  temperature  range  40  to  160  F  has  been  experimentally  determined 
(Ref.  42).  All  attempts  to  ut  ilize  the  estimation  techniques  available 
for  extrapolating  these  data  to  higher  temperatures  were  unsuccessful. 


Specific  Heat.  The  specific  heat  of  liquid  N^O^  from  20.5  to  64.8  F  was 
measured  by  Giauque  (Ref.  41 ) ,  Extrapolation  of  these  data  to  200  F 
was  made  by  means  of  Eq.  29. 

Cp  CO  0,680  =  0.0856  (29) 

Up  to  65  F,  the  results  obtained  are  identical  to  those  experimentally 
observed.  Above  200  F,  the  extrapolated  curve  shows  an  abrupt  change  of 
slope,  suggesting  that  the  equation  is  probably  not  applicable  at  these 
temperatures  (Fig.  24). 


Coefficient  of  Thermal  Expansion 

The  coefficients  of  thermal  expansion  for  N^O^  derived  from  the  experi¬ 
mental  densities  are  presented  in  Table  10. 
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Figure  24  .  Thermal  Conductivity  of  Liquid  NQ0 


,  btu/lb-f 


0.410 


Figure  25-  Specific  Heat  of  Liquid  N^O^ 
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TABLE  10 


COEFFICIENTS  OF  THERMAL  EXPANSION  FOR  NgO^ 


Temperature  Range,  F 

-ft  x  104/f 

32  to  70 

9.28 

70  to  130 

10.3 

130  to  190 

12.3 

190  to  250 

21.7 

250  to  310 

74.7 

SELECTION  OF  LIQUID  SIMULANTS 


The  liquid  simulants  for  the  selected  propellants  should  be  of  low  tox¬ 
icity,  nonflammable,  and  readily  available.  The  first  two  criteria  are 
of  importance  to  provide  maximum  safety  both  to  the  operators  conducting 
the  tests  and  to  the  equipment  employed.  Although  it  would  be  desirable 
to  use  the  simulant  for  repeated  tests,  it  should  be  available  in  rela¬ 
tively  large  quantities  and  require  a  minimum  cf  additional  processing. 

An  additional  factor  which  must  be  considered  is  the  condition  under 
which  the  propellant  will  be  simulated.  For  example,  if  the  structure 
and  propellant  to  be  tested  are  subjected  to  small  temperature  gradients 
or  slowly  occurring  temperature  changes,  even  distribution  of  heat  is 
allowed.  Therefore,  vapor  pressure  and  coefficient  of  thermal  expansion 
will  have  the  greatest  effect  on  the  structural  stresses.  Of  secondary 
importance  in  this  case  is  the  density.  For  the  simulant,  this  should  be 
no  less  than  that  of  the  propellant,  but  could  be  greater  to  increase 
the  g  effect  and  thus  check  the  safety  factor.  On  the  other  hand,  if 
the  propellants  are  subjected  to  large  temperature  gradients,  thermal 
conductivity,  viscosity,  and  specific  heat  become  more  important  because 
these  properties  will  determine  local  temperatures  and  temperature 
gradients . 

SIMULANTS  FOR  NO, 

2  k 

The  materials  which  are  recommended  for  use  as  simulants  for  N^O,  are 

2  4 

the  "Freons”  or  "Genetrons . "  The  compounds  of  particular  interest  are 
organic  compounds  which  may  contain  carbon,  hydrogen,  fluorine,  and 
chlorine.  The  wide  applicability  and  usefulness  of  these  compounds  are 
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due  to  the  unusual  combination  of  properties  found  in  them.  One  of  the 
most  important  qualities  of  these  fluorocarbons  is  the  low  level  of  tox¬ 
icity  (Ref.  44)  and  its  resultant  safety.  In  general,  the  Freons  exhibit 
excellent  thermal  and  chemical  stability.  None  are  flammable  or  explosive. 


Physical  Properties 

Four  of  the  Freons;  i.e.,  -11,  -21,  -112,  and  -113,  appear  to  be  service 
able  as  simulants  for  N^.  The  physical  properties  of  these  particular 
Freons  are  summarized  in  Table  11. 


Critical  Constants.  The  critical  constants  for  the  Freons  listed  in 
Table  12  are  measured  values  with  the  exception  of  the  critical  pressure 
for  Freon-112,  which  was  calculated. 


Vapor  Pressure.  The  vapor  pressure  for  Freon-11  (Fig.  26)  and  Freon-21 
(Fig-  27)  over  their  entire  liquid  range  was  calculated  by  use  of  E]q.  30 
and  31,  respectively. 

log  ~  =  l  (Tr)  -  0.285  0  (Tr)  (30) 

r 

log  jr  =  Jf  (Tr)  -  0.182  0  (Tr)  (31) 

r 

Experimental  values  for  the  vapor  pressure  of  Freon-112  (Fig,  28)  are 
available  in  the  temperature  range  50.6  to  194  F  (Ref.  45).  Extrapo¬ 
lation  of  the  data  to  the  critical  temperature  was  done  by  Eq.  32.  In 
those  cases  where  experimental  data  existed,  excellent  agreement  between 
the  calculated  and  experimental  values  was  observed. 

log  “  =  ft  (Tr)  +  0.232  0  (Tr)  (32) 

r 

The  vapor  pressure  data  for  Freon-113  is  shown  in  Fig.  29. 
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IY  OF  PHYSIC OCHEMICAL  PROPERTIES  OF  FKEONS  -11,  -21,  -112  and  -113 
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Density.  The  experimentally  determined  densities  of  Freon-11,  Freon-21, 
Freon-113  (Ref.  47),  and  Freon-112  (Ref.  45)  were  extrapolated  to  the 
critical  point  by  use  of  the  following  equations. 


Freon-11 : 

yP  (lb/ cu 

ft) 

=  760.7  co 

(Fig.  30) 

Fre on-21 : 

P  (lb/cu 

ft) 

=  719.0  CO 

(Fig.  31) 

Freon-112 : 

P  (lb/cu 

ft) 

=  796.6  CO 

(Fig.  32) 

Freon-113  * 

P  (lb/cu 

ft) 

=  1.6212 

2.172  X 

-  3-3  x  10-6t2 


(Fig.  33) 

where  t  is  in  centigrade  degrees 


In  all  cases,  the  deviations  between  experimental  and  calculated  values 
were  less  than  0.5  percent. 


Viscosity 

A  search  of  the  literature  did  not  yield  any  significant  amount  of  data 
on  the  viscosity  of  the  Freons.  Only  for  Freon-112,  were  viscosity  data 
found  over  the  temperature  range  of  interest  (Ref.  45).  Therefore,  for 
Freon-11,  -21,  and  -113,  the  method  of  Uyehara  and  Watson  was  used  for 
the  calculation  of  viscosity  as  a  function  of  temperature.  The  following 
equations  were  employed: 


Freon  - 

11: 

V 

(lb/ft-sec) 

=  2.78  x 

101  ’’r 

(Fig.  34) 

Freon  - 

21: 

n 

(lb/ f t-sec ) 

=  2.58  x 

—5 

10  J  -n 
'r 

(Fig.  35) 

Freon  - 

112: 

V 

(lb/ft-sec) 

1 

0 

00 

X 

IQ"5  r\ 

(Fig.  36) 

Freon  - 

113: 

V 

(lb/ft-sec) 

=  3-85  x 

O 

1 

(Fig.  37) 

Thermal  Conductivity.  Experimental  data  have  been  reported  for  the  ther¬ 
mal  conductivity  of  Freon  -11,  -21  and  -113  in  the  temperature  range  32 
to  167  F  (Ref.  48).  The  experimental  data  were  extrapolated  to  the  critical 
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DENSITY,  LB/CU-FT 


TEMPERATURE, F 


Figure  30.  Orthobaric  Density  of  Freon-11 
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DENSITY,  LB/CU-FT 


TEMPERATURE,  F 


Figure  31.  Orthobaric  Density  of  Freon-21 
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Figure  32.  Orthobaric  Density  of  Freon-112 
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Figure  34.  Viscosity  of  Freon-11 
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TIC  VISCOSITY, CENTISTOKES 


VISCOSITY,  LB/FT-SEC 


temperature,  f 

Figure  35.  Viscosity  of  Freon-21 
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KINEMATIC  VISCOSITY,  CENT I STOKES 


Figure  37.  Viscosity  of  Freon-113 
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KINEMATIC  VISCOSITY,  CENTISTOKES 


point  by  use  of  the  equations  given  in  Ref. 48  .  For  Freon  -112,  Weber's 
equation  was  used  to  estimate  thermal  conductivity. 


The  general  equation  given  for  Freon  -11,  -21,  and  -113  for  the  tempera¬ 
ture  range  32  to  167  F  is: 


K, 


*b 


1  +  a  (t  -  t  ) 


(33) 


where 


“b 

a 


thermal  conductivity  in  Btu/hr-ft-F  at  the  desired 
temperature  (F) 

thermal  conductivity  in  Btu/hr-ft-F  at  32  F 

the  temperature  coefficient  of  thermal  conductivity 


The  following  values  of  K_  and  Of  for  the  various  Freons  were  used. 


Freon 

*b 

a  X  103 

11 

0.0680 

-1.928 

21 

0.0770 

-1.73 

113 

0.0576 

-1.745 

The  thermal  conductivities  of  the  four  Freons  are  given  in  Fig. 38  through 
41. 


Heat  Capacity.  The  experimentally  determined  heat  capacities  of  the 
Freons  (Ref.  46,  47,  and  49)  were  extrapolated  by  the  following  equations 
to  a  temperature  at  which  rapid  change  in  the  slope  was  observed. 
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Figure  38.  Thermal  Conductivity  of  Freon-11 
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Figure  39.  Thermal  Conductivity  of  Freon-21 
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Figure  41.  Thermal  Conductivity  of  Freon-113 
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Freon- 11  : 
Freon-21  : 
Fre on-1 12 : 
Fre on-1 13 : 


C 

P 

C 

P 

C 

P 

C 

P 


0.0568  (Fig.  42) 
0.1066  (Fig.  43) 
0.0406  (Fig.  44) 
0.0245  (Fig.  45) 


Evaluation  of  the  data  collected  led  to  the  conclusion  that  Freon  -21 
(Table  12)  would  serve  best  as  a  simulant  for  throughout  the  liquid 

temperature  range  of  the  propellant.  This  appears  to  be  the  case  even 
though,  up  to  approximately  150  F,  the  vapor  pressure  of  Freon  -21  is 


slightly  higher  than  that  of  N^O^.  Above  this  temperature,  the  vapor 


pressure  of  the  Freon  becomes  less  than  that  of 

ties  of  Freon  -21  resemble  more  nearly  those  of  NO, 

J  2  4 


Freons . 


All  other  proper- 
than  do  the  other 


The  second  Freon  chosen  as  a  simulant  for  N^O,  is  Freon  -11.  In  this 

2  4 

case,  the  heat  capacity  of  Freon  -11  over  the  entire  liquid  temperature 
range  of  is  considerably  less  than  that  of  the  propellant. 

Because  their  densities  are  somewhat  higher  than  throughout  the 

liquid  temperature,  Freons  -112  and  -115  could  be  used  to  test  the 
effect  of  g  loading.  However,  for  more  complete  testing  utilizing  only 
one  simulant,  the  properties  of  Freons  -112  and  -113  are  not  as  close  to 
those  of  as  are  those  of  Freon  -11  and  Freon  -21.  Finally,  the  high 

melting  point  of  Freon  -112  (79  F)  precludes  its  use  in  the  lower  tem¬ 
perature  range. 
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A  COMPARISON  OF  THE  PHYSICAL  PROPERTIES 
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'ABLE  12 
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10.3  9-9  10.9  8.8 

74.7  20.6  29.6  16.6 


0.30 


Heat  of  Freon-11  Figure  43.  Specific  Heat  of  Freon-21 


0.30 
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Figure  44.  Specific  Heat  of  Freon-112  Figure  45.  Specific  Heat  of  Freon-113 


SIMULANTS  FOR  MON (METHYLH YDRA Z I NE  AND 
TITAN  II  FUEL  BLEND 


The  similarity  of  both  the  chemical  composition  and  the  physical  prop¬ 
erties  of  monomethylhydrazine  and  the  Titan  II  fuel  blend  suggested  that 
a  simulant  suitable  for  one  might  also  serve  for  the  other.  Because  of 
the  polar  character  of  the  solvents,  simulants  which  were  also  polar  were 
considered  and  investigated.  Among  the  compounds  and  mixtures  which  were 
considered  were  water,  cellosolves,  carbitols,  and  water-alcohol  mixtures 
In  the  following  discussions,  only  those  materials  which  appear  to  be  the 
best  simulants  are  described. 


Water 


Water,  because  of  its  nonhazardous  nature  and  availability,  would  be  an 
ideal  simulant.  In  general,  it  appears  to  be  a  suitable  simulant  for 
the  Titan  II  fuel  blend  through  nearly  all  of  its  entire  liquid  range, 
and  would  similarly  serve  as  a  simulant  for  monomethylhydrazine  through 
most  of  its  liquid  range. 


Physical  Properties.  A  summary  of  the  physical  properties  of  water  is 
given  in  Table  13  along  with  those  of  the  other  suggested  simulants. 

The  graphs  showing  the  variance  of  a  property  with  temperature  were  pre¬ 
pared  from  experimental  data.  Equations  that  approximate  the  experi¬ 
mental  data  to  the  degree  of  accuracy  indicated  are  presented. 
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A  COMPARISON  OF  THE  PHYSICAL  PROPERTIES  OF  MON OMETHYLHYDRAZ INE 
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TABLE  13 
(Continued ) 
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Vapor  Pressure.  The  vapor  pressure  of  water  has  been  determined 
over  the  entire  temperature  range  of  32  to  705.  6  F  (Bef.  50).  The 
results  (Fig.  46)  are  satisfactorily  represented  (average  deviation  3*5 
percent  from  1  to  3200  psia)  by  the  equation: 

log  =  0  (Tr)  +  0.411  (Tr)  (34) 

r 

Density.  The  density  of  water  as  a  function  of  temperature  (Fig.  47) 
over  its  entire  liquid  range  (Bef.  50)  is  correlated  by  Eq.35  with  an 
average  deviation  of  3*1  percent. 

p  (lb/cu  ft)  =  474.4  «  (35) 


Viscosity.  The  viscosity  of  water  (Fig.  48)  expressed  as  a  function 
of  temperature  (Ref.  50  and  5l)  is  governed  by  Eq.36  (Bef.  5l). 


log  r\ 

=  A  +  BT  +  ^ 

where 

V  = 

millipoise 

A  = 

-2.5608 

B  - 

-1.627  x  10"5 

C  = 

120.63 

D  - 

-176.98 

T  - 

absolute  temperature 

(36) 


Thermal  Conductivity.  The  thermal  conductivity  of  water  (Fig. 49  ) 
is  known  up  to  690  F  (Ref.  52).  Smith's  equation  (Bef.  ll)  approximates 
the  data  to  an  average  deviation  of  13  percent. 
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Figure  47.  Orthobaric  Density  of  Water 
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Figure  48.  Viscosity  of  Water 
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KEJEMATIC  VISCOSITY  (CEIITISTOKES) 
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Figure  49.  Thermal  Conductivity  of  Water 


96 


Heat  Capacity.  The  heat  capacity  data  for  water  (Fig.  50  )  is  com¬ 
plete  up  to  690  F  (Ref.  52).  Equation  37  >  derived  from  the  Chow-Bright 
correlation,  represents  the  data  within  an  average  deviation  of  less 
than  1  percent. 

C  CO0-06  =  0.886  (37) 

P 

Coefficient  of  Thermal  Expansion.  Coefficients  of  thermal  expan¬ 
sion  for  water  are  given  in  Table  14. 


TABLE  14 

COEFFICIENTS  OF  THERMAL  EXPANSION  OF  WATER 


Temperature  Range,  F  x  104/F 
32  to  212  2.39 

212  to  356  5-67 
356  to  500  9.13 
500  to  644  19.9 


A  comparison  of  the  physical  properties  of  MMH  and  the  Titan  II  fuel 
blend  with  those  of  water  at  temperatures  of  32  to  approximately  600  F 
indicates  the  following  (see  Table  13). 

Vapor  Pressure.  The  vapor  pressure  of  water  is  lower  than  that  of 
the  50/50  fuel  blend  throughout  the  latter1 s  liquid  range,  and  is  also 
lower  than  that  of  MMH  up  to  a  temperature  of  550  F. 
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TEMPERATURE,  F 

Figure  50.  Heat  Capacity  of  Water 
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Density.  The  density  of  water  is  greater  than  that  of  the  50/50 
fuel  blend,  as  well  as  that  of  MMH  over  the  liquid  range  of  water.  This 
provides  a  means  of  checking  the  g  loading. 


Viscosity.  The  absolute  and  kinematic  viscosities  of  water  are 
slightly  greater  than  those  of  the  two  propellants  in  the  temperature 
range  32  to  100  F.  Above  100  F,  the  viscosities  are  only  slightly  less 
than  that  of  the  two  propellants. 


Thermal  Conductivity  and  Heat  Capacity.  The  simulation  of  the 
thermal  conductivity  and  heat  capacity  of  MMH  and  the  50/50  fuel  blend 
by  water  is  not  completely  satisfactory.  In  the  case  of  the  heat  capac¬ 
ity,  the  data  for  water  are  approximately  50  percent  higher  than  that  of 
the  two  propellants.  The  positive  deviation  of  the  thermal  conductivity 
values  is  even  greater.  The  limited  heat  capacity  data  available  for  the 
propellants  at  the  lower  temperatures  indicates  that  the  deviation  is 
real  and  not  solely  the  consequence  of  the  extrapolation  methods  used. 

In  view  of  the  measured  thermal  conductivity  values  for  the  50/50  fuel 
blend,  it  is  most  likely  that  the  deviations  observed  are  also  real  but 
the  actual  extent  is  not  accurately  known. 


Carbitols  and  Cell os ol ve s 


The  physical  properties  of  various  carbitols  and  cellosolves  were  investi¬ 
gated  as  possible  simulants  for  MMH  and  the  Titan  II  fuel  blend.  The 
materials  which  approximate  all  of  the  physical  properties  of  the  propel¬ 
lants  best  are  diethyl-  and  dibutylcarbitol  and  butylcellosol ve .  In  the 
case  of  these  simulants,  it  was  found  that  only  the  calculated  thermal 
conductivity  differed  significantly. 


99 


Physicochemical  Properties.  A  summary  of  the  pertinent  physical  prop¬ 
erties  of  the  selected  carbitols  and  cellosolves  is  given  in  Table  15- 

Critical  Constants.  No  experimental  data  on  critical  constants  of 
the  selected  carbitols  and  cellosolves  were  found.  The  critical  pres¬ 
sure  and  temperature  were  therefore  calculated  by  means  of  Lydersen's 
method  (Ref.  1  ). 

Vapor  Pressure.  The  vapor  pressure  of  diethyl-  and  dibutylcarbitol 
was  calculated  from  the  boiling  points  of  the  material.  The  vapor  pres¬ 
sure  of  butylcellosolve  has  been  experimentally  determined  over  a  short 
temperature  range  (Ref.  53)*  Extrapolation  of  the  vapor  pressure  was 
made  using  the  equation  given  in  Ref.  53-  The  vapor  pressures  are  plot¬ 
ted  as  functions  of  the  temperature  in  Fig. 53  through  56. 

Diethylcarbitol :  log  — -  =  0  (T^)  +  2.045  0  (T^)  (38) 

r 

Dibutylcarbitol:  log  =  jA  (T^)  +  3-096  if)  (T^)  (39) 

r 

Density.  The  experimental  densities  of  diethylcarbitol  (Ref. 54), 
dibutylcarbitol  (Ref.  55),  and  butylcellosolve  (Ref.  56)  were  extrapo¬ 
lated  by  means  of  Eq .  40,  41,  and  42,  respectively.  The  data  are  pre¬ 
sented  in  graphical  form  in  Fig. 54  through  56. 

p  (lb/cu  ft)  =  417-3  ^ 

p  (lb/cu  ft)  =  398.9  td 

p  (lb/cu  ft)  -  415.8  00 


100 


(40) 

(41) 

(42) 


SUMMARY  OF  PHYSICOCHEMICAL  PROPERTIES  OF  WATER,  AM)  SELECTED 
CARBITOLS  AND  CELLOS OLVES 
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Figure  51.  Vapor  Pressure  of  Diethyl  Carbitol 
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Figure  54.  Orthobar ic  Density  of  Diethyl  Carbitol 
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Viscosity.  The  absolute  viscosities  of  the  two  carbitols  were  esti¬ 
mated  (Ref.  57)  and  extrapolated  by  the  method  of  Uyehara  and  Watson.  The 
experimental  viscosity  data  of  butylcellosolve  (Ref .58  )  was  extrapolated 
via  the  Uyehara  and  Watson  equation.  The  following  equations  for  estima¬ 
tion  of  the  viscosities  were  derived: 

Diethylcarbitol :  T?  (lb/ft-sec)  =  2.6l  x  10  ^  (Fig.  57)  (43) 

Dibutylcarbitol :  V  (lb/ft-sec)  =  1.68  x  10  ^  (Fig.  58)  (44) 

Butylcellosolve:  r?  (lb/ft-sec)  =  4.89  x  10"5  r?r  (Fig.  59)  (45) 

The  extrapolated  data  for  butylcellosolve  was  checked  at  100,  210,  and 
250  F.  The  measured  kinematic  viscosity  of  butylcellosolve  at  100  F  was 
2.37  centistokes  (calculated,  2.24  centistokes) ;  at  210  F,  0.84  centi- 
stokes  (calculated,  1.23  centistokes);  and  at  250  F,  O.63  centistokes 
(calculated,  1.06  centistokes). 

These  measurements  seem  to  indicate  that  for  the  cellosolves  (compounds 
which  can  undergo  hydrogen  bonding  readily)  the  estimation  technique  is 
incapable  of  accounting  for  the  change  in  the  degree  of  aggregation  with 
increases  in  temperature.  It  is  therefore  expected  that  the  calculated 
viscosities  of  butylcellosolve  at  the  higher  temperatures  may  be  in  error 
by  perhaps  as  much  as  100  percent.  However,  for  the  carbitols,  because 
hydrogen  bonding  cannot  occur,  the  estimation  of  the  viscosities  should 
be  significantly  better. 


Heat  Capacity.  The  heat  capacity  of  the  carbitols  was  estimated 
by  the  method  of  Johnson  and  Huang  (Ref. 57  ).  This  method  has  been 
found  to  give  heat  capacities  at  68  F  with  an  error  of  ±2  to  4  percent. 
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Figure  57.  Viscosity  of  Diethyl  Carbitol 
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Figure  58. Viscosity  of  Dibutyl  Carbitol 
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The  estimated  heat  capacity  of  the  compound  was  extrapolated  by  the 
method  of  Chow  and  Bright,  using  the  value  of  2.8  for  the  constant  "a," 
and  determining  the  constant  Mb"  for  each  compound  from  the  single  heat 
capacity  value  obtained  from  the  method  of  Johnson  and  Huang. 

The  experimental  heat  capacity  of  butylcellosolve  was  extrapolated  using 
the  method  of  Chow  and  Bright.  The  following  equations  relate  the  heat 
capacity  data,  plotted  in  Fig.  60  through  62,  with  temperature. 


Diethylcarbitol : 

c  CO2-8  = 
p 

0.00188  (Fig.  60) 

(46) 

Dibutylcarbitol : 

C  O)2-8  = 

p 

0.00194  (Fig.  61) 

(47) 

Butylcellosolve : 

c  ft)1-4*. 

p 

0.03178  (Fig.  62) 

(48) 

Thermal  Conductivity.  Thermal  conductivity  data  for  the  three 
simulants  are  available  in  the  temperature  range  86  to  194  F  (Ref.  58). 
Extrapolation  of  the  experimental  data  was  performed  by  means  of  the 
equations  given  in  Ref.  58  for  this  temperature  range.  The  extrapolated 
thermal  conductivities  would  be  expected  to  be  only  crude  estimates  as 
the  temperature  approaches  the  critical  temperature.  The  thermal  con¬ 
ductivity  data  are  given  in  Fig.  63  through  65. 


Evaluating  all  the  data  obtained  for  the  four  candidate  simulants,  it 
appears  that  the  experimental  and  extrapolated  values  of  the  heat  capac¬ 
ities  and  thermal  conductivities  of  diethylcarbitol ,  dibutylcarbitol  and 
butylcellosolve  agree  more  closely  with  those  of  MMH  and  the  Titan  II 
fuel  blend  than  do  those  of  water  (Table  15).  The  other  properties  of 
the  propellants  such  as  viscosity  and  density  are  simulated  very  well  by 
all  three  materials.  In  addition,  the  vapor  pressures  of  the  three 
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Figure  60.  Specific  Heat  of  Diethyl  Carbitol 
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Figure  6l .  Specific  Heat  of  Dibutyl  Carbitol 
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Figure  62. Specific  Heat  of  Butyl  Cellosolve 
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Figure  63.  Thermal  Conductivity  Diethyl  Carhitol 
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Figure  64.  Thermal  Conductivity  of  Dibutyl  Carhitol 
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Figure  65. Thermal  Conductivity  of  Butyl  Cellosolve 
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materials  are  significantly  lower  than  those  of  the  propellants.  Of  the 
three  compounds,  the  best  simulant  is  dibutylcarbitol .  The  factors  which 
governed  this  choice  are:  the  low  melting  point,  the  relatively  high 
flash  point  of  245  F,  and  the  over-all  better  agreement  of  its  physical 
properties  with  those  of  the  two  propellants.  The  second  and  third 
choices  are  diethylcarbitol  and  butylcellosolve ,  respectively. 


SIMULANTS  FOR  RP-1  AND  JP-4 

Three  materials:  dibutylcarbitol,  red  oil  (commercial  oleic  acid),  and 
olive  oil  were  evaluated  as  simulants  for  RP-1  and  JP-4.  Butylcello¬ 
solve  was  another  candidate  but  has  been  ruled  out  because  of  the  lower 
flash  point  (l60  F).  The  properties  of  dibutylcarbitol  have  been  reported 
previously  (Table  13). 


Red  Oil  and  Olive  Oil 


Red  oil  and  olive  oil  have  high  flash  points  of  372  and  437  F,  respect¬ 
ively  which  makes  these  materials  attractive  as  simulants  for  RP-1  at 
temperatures  which  might  exceed  the  flash  point  of  dibutylcarbitol. 

A  summary  of  the  physicochemical  properties  of  the  two  oils  is  given 
in  Table  16.  The  equations  which  were  used  in  extrapolating  or  esti¬ 
mating  the  physical  properties  of  the  oils  are  given  in  the  following 
sections.  Where  experimental  data  were  found  for  a  specific  physical 
property,  these  are  reported  also.  The  figure  number  which  corresponds 
to  the  graph  relating  physical  property  with  temperature  is  also 
indicated . 
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TABLE  16 
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Vapor  Pressure. 


1 

=  ji  (Tr)  +  16.22  0  (Tr) 

(Ref.  60) 

Olive  oil:  log 

1 

P 

(Fig.  66) 

Red  oil:  log 

r 

_i_ 

p 

r 

=  ji  (Tr)  +  4.55  0  (Tr) 

(Ref.  6l) 

(Fig.  67) 

Density. 

Olive  oil:  p  (lb/cu  ft)  =  391*6  CO 
Red  oil:  p  (lb/cu  ft)  =  389-7  <0 


(Ref.  62) 

(Fig.  68) 

(Ref.  63  and  64) 
(Fig.  69; 


Viscosity. 

Olive  oil: 
Red  oil: 


TJ  (lb/ft-sec)  =  0.1418  x  10"4  rjr 
r\  (lb/ft-sec)  =  0.569  x  10~4 


Ref.  55) 
Fig.  70) 

Ref.  64) 
Fig.  7l) 


The  values  calculated  from  the  above  two  equations  are  slightly  higher  than 
those  from  the  graph.  The  curves  shown  in  Fig.  70  and  71  were  drawn  to 
reflect  the  general  curvature  of  the  experimental  viscosity  data  at  the 
lower  temperatures. 

Heat  Capacity. 


Olive  oil: 

2  8  -4 

C  <J)*'  =  2.18  x  10 

P 

(Ref. 

(Fig. 

55) 

72) 

Red  oil: 

3  3  -4 

C  td  =  4.90  x  10 

(Ref. 

64) 

P 

(Fig. 

73) 

Thermal  Conductivity. 

Olive  oil: 

Red  oil: 


K(Btu/hr-ft-F>  =  9.8  x  10 
K(Btu/hr-ft-F)  =  1.089  Cp  p 


r 

L 


1-3.5x10 


(o) 

M 


.  (Ref.  65, 
t|  66,  and  67) 
(Fig.  74) 

Ref.  64) 

Fig.  75) 
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PRESSURE,  PSIA 


10 


TEMPERATURE,  F 


Figure  66.  Vapor  Pressure  of  Olive  Oil 
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KINEMATIC  VISCOSITY,  CENT I STOKES 


PRESSURE 


TEMPERATURE,  F 

Figure  67.  Vapor  Pressure  of  Red  Oil 
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p ,  lb/cu-ft 


60 


TEMPERATURE,  F 

Figure  68.  Density  of  Olive  Oil 
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1C  VISCOSITY,  CENT I  STOKES 


Figure  71.  Viscosity  of  Red  Oil 
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KINEMATIC  VISCOSITY,  CENT  1ST OKES 


,  btu/lb-f 


TEMPERATURE,  F 

Figure  72. Heat  Capacity  of  Olive  Oil 
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Figure  73. Heat  Capacity  of  Red  Oil 
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Figure  75.  Thermal  Conductivity  of  Red  Oil 
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An  evaluation  of  all  the  data  compiled  (Table  17  and  17A)  shows  that  the 
physical  properties  of  RP-1  are  simulated  quite  well  by  dibutylcarbitol. 
The  one  point  which  is  not  entirely  satisfactory  is  the  flash  point  of 
dibutylcarbitol,  which  is  only  about  100  degrees  higher  than  that  of  RP-1. 

Estimation  of  the  physical  properties  of  red  oil  and  olive  oil  at  various 
temperatures  indicates  that  the  main  drawback  to  the  use  of  these  materi¬ 
als  as  simulants  for  HP-1  is  their  significantly  greater  viscosities. 

All  other  properties  of  these  materials  are  in  good  agreement  with  those 
of  RP-1  at  the  elevated  temperatures.  The  relatively  high  melting  point 
of  63  F  for  red  oil  limits  its  usefulness  as  a  simulant  for  RP-1  to  tem¬ 
peratures  above  which  olive  oil  decomposes  (approximately  330  F). 

The  high  viscosity  of  the  vegetable  oils  may  serve  as  a  more  severe  test 
in  temperature  gradients  at  the  wall  because  heat  cannot  be  carried  away 
as  rapidly  as  for  the  propellant  itself.  It  has  been  reported  (Ref.  7  ) 
that  olive  oil,  when  heated  for  2-3/4  hours  at  6 1 5  F,  is  50.4-percent 
decomposed.  At  this  temperature,  the  oil  undergoes  cracking  and  poly¬ 
merization.  This  leads  ultimately  to  the  formation  of  cross-linked  solid 
resins . 

In  the  case  of  JP-4,  the  temperature  range  of  interest  to  the  Air  Force 
was  -60  to  +230  F.  For  the  major  portion  of  this  range  (up  to  approxi¬ 
mately  125  to  150  F),  dibutylcarbitol  will  be  satisfactory.  Above  this, 
the  vegetable  oils  are  recommended,  with  red  oil  the  material  of  choice. 
Selection  of  this  oil  is  based  on  the  low  decomposition  temperature  of 
olive  oi  1  (330  F).  Red  oil  can  probably  be  employed  from  above  its  melt¬ 
ing  point  (63  F)  up  to  the  pseudocritical  temperature  of  JP-4. 
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TABLE  17 


A  COMPARISON  OF  THE  PHYSICAL  PROPERTIES  OF 
RP-1  WITH  ITS  SIMULANTS 


T,  F 

RP-l 

Dibuty  1  carb  i  tal 

Olive  Oil 

Red  Oil 

U  1  ■  1  1  II  III  1  Wl  !!■■■■—— 1 

100 

0.06 

Nil 

Nil 

Nil 

4.5 

0.4 

Nil 

0.02 

■99 

45 

17 

Nil 

1.0 

Viscosity,  X  x  10^  lb/ft-sec 

100 

90 

70 

2300 

23 

18 

140 

130 

10 

8 

35 

37 

Heat  Capacity,  Btu/lb-F 

100 

0.48 

0.52 

0.50 

0.46 

0.58 

0.70 

0.62 

0.64 

0.68 

1.09 

0.79 

0.90 

Density,  lb/cu  ft 

100 

49.4 

54.3 

56.3 

55.0 

44.4 

48.8 

52.0 

50.5 

■9 

38.3 

41.8 

47.5 

45.3 

Thermal  Conductivity,  Btu/hr-ft-F 

mm 

1 

0.083 

0.097 

1 

1 

0.095 

0.077 

0.072 

0.093 

0.078 
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TABLE  17A 


A  COMPARISON  OF  JP-4  WITH  ITS  SIMUIANTS 
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SIMULATION  OF  SOLID  PROPELLANTS 


BACKGROUND 

Unlike  liquid  propellants  and  ingredients  which  exist  separately  in  a 
missile,  a  solid  propellant  and  its  components  are  intimately  mixed  and 
in  contact  with  each  other  in  the  vehicle.  Thus,  successful  simulation 
is  not  merely  a  matter  of  duplication  of  the  properties  of  the  individual 
ingredient,  but  rather  a  very  exact  matching  of  the  properties  of  the  pro¬ 
pellant  formulation.  This  duplication  of  properties  should  be  attained 
without  any  changes  from  the  compositional  content  of  the  live  propellant. 
As  with  the  liquids,  satisfactory  simulation  was  based  on:  (l)  safety  in 
use  (as  reflected  by  little  or  no  fire  or  toxicity  hazard),  and  (2)  dupli¬ 
cation  of  the  mechanical  and  thermal  properties  of  the  live  propellant. 

There  has  been  a  long  history  of  attempts  at  simulation  of  solid  propel¬ 
lant  compositions.  However,  most  of  this  work  has  been  done  either  to 
check  out  processing  equipment  with  inert  compositions,  to  study  and 
develop  new  and  improved  processing  procedures,  or  to  train  personnel. 

In  most  instances  of  composite  propellant  simulation,  the  rheological 
properties  of  the  simulated  formulations  have  not  been  sufficiently  close 
to  those  of  the  live  compositions  for  satisfactory  translation  of  results 
obtained  with  the  former  types  into  improved  procedures  for  processing 
the  latter. 

In  the  past,  simulants  for  ammonium  perchlorate  have  included  glass 
spheres,  sodium  chloride  or  sulfate,  and  ammonium  chloride  or  sulfate. 

For  a  large  portion  of  the  work  done  at  the  Naval  Ordnance  Test  Station 
on  the  screw  extrusion  of  double  base  propellant,  an  inert  composition 
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based  on  cellulose  acetate  and  a  phthalate  or  azelate  plasticizer  was 
employed  to  simulate  rheological  properties.  A  considerable  amount  of 
information  on  the  effect  of  oxidizer  particle  size  and  size  distribution 
has  been  obtained  in  simulated  systems  using  glass  beads.  However,  the 
data  obtained  did  not  agree  quantitatively  with  that  from  live  formulations. 

Similarly,  the  limited  amount  of  work  that  has  been  done  on  the  mechanical 
properties  of  propellants  via  simulation  studies  has  not  been  particularly 
successful  in  duplicating  the  properties  of  the  real  composition  because 
of  the  following  factors: 

1.  The  polymer-oxidizer  bond  strength 

2.  The  effect  of  ions  of  dissolved  oxidizer  on  the  curing  of  the 
polymer 

3.  The  solubility  interaction  of  the  polymer,  oxidizer,  and 
plasticizer 


Polymer-Oxidizer  Bond  Strength 

An  examination  of  the  actual  nature  of  a  composite  propellant  reveals 
that  very  little  of  the  polymer  can  be  more  than  one  micron  away  from  a 
solid  surface.  This  stems  from  optimum  packing  consideration.  For  poly¬ 
mer  chains  of  a  molecular  weight  of  2000,  it  is  likely  that  10  percent  of 
these  units  actually  touch  a  surface  and  probably  all  of  the  units  are 
within  the  electrostatic  field  of  the  surface  Van  der  Waal  forces.  In 
many  composite  propellants,  the  extent  of  these  surface  interactions 
determines  the  mechanical  properties  of  the  cured  propellants  to  a  larger 
degree  than  do  the  binder* s  intrinsic  mechanical  properties.  The  higher 
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the  polarity  of  the  binder,  the  greater  the  surface  interaction  for  the 
same  oxidizer.  For  the  same  binder,  the  higher  the  ionic  strength  of  the 
oxidizer,  the  greater  the  interaction.  Thus,  it  may  be  seen  that  highly 
ionic  alkali  metal  salts  should  not  be  substituted  for  ammonium  salts,  or 
vice  versa.  In  addition,  it  is  not  unexpected  that  glass  beads  would 
serve  as  poor  simulants  for  ionic  substances  such  as  the  oxidizer  salts. 


Solubility  of  Ingredients 

Ammonium  perchlorate  dissolves  in  a  very  large  number  of  materials  which 
are  polar  to  any  extent.  The  trace  solubility  per  se  is  of  little  import 
by  itself,  but  the  dissolved  ammonium  ion  has  a  very  marked  effect  on  the 
curing  reactions  of  polyurethanes  and  butadiene  copolymers.  These  curing 
reactions  occur  between  an  active  hydrogen  and  some  other  active  group. 

In  the  former  case,  urethane  formation  is  affected  by  the  ammonium  ion; 
for  the  latter,  the  opening  of  the  oxirane  or  azirane  ring  of  the  curing 
agent  is  also  influenced  by  the  dissolved  ingredient.  The  ammonium  ion 
can  act  either  as  a  competitive  proton  donor  or  as  a  suppressor  to  the 
activity  level  of  the  active  hydrogen  in  the  monomer  unit.  Even  differ¬ 
ent  ammonium  salts  have  different  effects  on  the  rates  and  degree  of  cure. 
The  apparently  nearly  analogous  hydrazinium  perchlorate  has  a  markedly 
different  effect  on  the  curing  reactions  than  the  ammonium  compound. 

Thus,  it  is  not  surprising  that  so  different  a  material  as  a  metal  salt 
is  a  very  poor  simulant  for  an  ammonium  salt. 

Polymer-Plasticizer-Oxidizer  Interactions 

In  the  case  of  plasticized  binder  systems  such  as  double-base  compositions 
or  cases  where  diluents  such  as  low-molecular-weight  hydrocarbons  or  nitro 
containing  ingredients  are  added  to  composite  systems,  a  three-way  set  of 
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interactions  must  be  considered.  The  binder-oxidizer  bond  strength  dis¬ 
cussed  above  is  further  complicated  by  a  third  component  (the  plasticizer 
or  diluent)  whose  presence  changes  the  intrinsic  binder  properties  mark¬ 
edly.  In  addition,  the  binder-oxidizer  bond  strength  will  vary  widely 
with  the  degree  of  plasticizer  accumulation  at  the  oxidizer  surfaces. 
Simulation  of  this  type  of  interaction  requires  a  good  understanding  of 
solubility  parameters  and  partition  phenomena. 

From  the  discussion  above,  it  can  be  seen  that  not  only  are  there  marked 
effects  of  surface  interactions,  but  that  the  cured  polymer  formed  in 
the  presence  of  included  solids  has  a  decidedly  different  structure  from 
that  formed  in  the  absence  of  filler.  In  fact,  the  polymer  structure 
probably  differs  from  one  filler  to  another.  Thus,  it  is  not  surprising 
that  attempts  to  predict  propellant  behavior  from  binder  properties  have 
been  unsuccessful . 


ROCKETDYNE  APPROACH  TO  THE  PROBLEM 


The  solid  propellants  to  be  simulated  were  aluminized  formulations  con¬ 
taining  ammonium  perchlorate  and  a  polybutadiene  resin. 

It  was  concluded  that  the  simulation  of  the  ammonium  perchlorate  must 
almost  of  necessity  be  done  by  an  ammonium  salt.  The  use  of  an  ammonium 
salt  of  a  strong  acid  such  as  ammonium  sulfate  would  provide  a  stable 
material  of  similar  properties  to  ammonium  perchlorate.  The  differences 
from  one  to  the  other  were  assumed  to  be  due  to  the  acidities  of  the  dis¬ 
solved  salts  in  the  binder  and  the  salt-binder  interfacial  bond  strength. 
This  latter  effect  is  believed  to  be  far  more  important  than  the  crystal 
structure,  particularly  where  the  simulant  has  approximately  the  same 
dimensions  and  shape  as  ammonium  perchlorate. 


138 


Total  replacement  of  perchlorate  by  sulfate  was  ruled  out  because  previ¬ 
ous  attempts  had  not  resulted  in  satisfactory  duplication  of  properties. 

The  approach  employed  was  the  substitution  of  ammonium  sulfate  for  the 
perchlorate  in  the  coarse-size  oxidizer,  while  retaining  the  perchlorate 
for  the  fine  size.  The  higher  specific  surface  of  the  fine  size  results 
in  this  fraction  having  from  two  to  ten  times  the  surface  of  the  coarse 
fraction,  even  though  the  fine  size  makes  only  a  small  weight  contribu¬ 
tion.  Therefore,  the  fine  size  would  be  the  major  source  of  the  dis¬ 
solved  ammonium  salts  and  the  resulting  surface  interaction.  The  coarse 
fraction,  however,  represents  60  to  80  percent  of  the  mass  of  the  oxi¬ 
dizer.  Therefore,  the  substitution  of  ammonium  sulfate  for  the  coarse 
perchlorate  reduces  the  oxidizer  to  a  fraction  of  its  original  level. 

In  the  presence  of  the  mixed  oxidizer  simulant,  the  binder  reacts  as  if 
it  is  in  contact  solely  with  ammonium  perchlorate. 

Using  this  approach,  the  simulation  of  RDS  501,  an  aluminized  propel¬ 
lant  containing  ammonium  perchlorate,  and  a  butadiene  resin  (Phillips 
Butarez  CTL)  was  evaluated  in  regard  to  safety  and  the  duplication  of 
mechanical  and  thermal  properties.  The  composition  of  RDS  501  is 
described  in  Appendix  A. 


Safety  Tests 

Several  years  ago  it  was  observed  at  Rocketdyne  that  a  propellant  which 
contains  only  ammonium  sulfate  in  the  coarse-oxidizer  phase  would  not 
ignite  under  ambient  conditions  even  when  a  blow  torch  flame  was  applied. 
At  the  start  of  this  program,  compositions  similar  to  RDS  501,  but  in 
which  50,  75,  and  100  percent  of  the  coarse  perchlorate  was  replaced  by 
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sulfate,  were  prepared  and  cured.  Strands,  0.25  inch  in  diameter  and 
3  to  4  inches  in  length,  were  cut  from  each  batch.  Ignition  wires 
embedded  in  live  propellant  were  bonded  to  each  end  of  the  individual 
strands.  Table  18  reviews  the  results  of  these  "hot-wire11  ignition 
tests.  The  only  composition  which  would  not  ignite  at  elevated  pres¬ 
sures  was  the  one  in  which  the  total  replacement  of  the  coarse  per¬ 
chlorate  by  sulfate  had  been  made.  Figure  ?6  is  a  photograph  of  sam¬ 
ples  of  such  a  composition  after  ignition  testing.  Only  minor  charring 
occurred  even  at  pressures  as  high  as  1440  psi.  In  a  somewhat  more 
severe  test,  a  strand  was  shredded  (to  provide  greater  burning  surface). 
It  still  did  not  ignite  under  ambient  conditions  (Fig.  77  ,  far  right). 

As  a  final  demonstration  of  the  safety  of  the  simulant,  a  gas  generator 
igniter  was  employed  rather  than  the  hot  wire  used  in  the  previous  tests. 
Again,  ignition  did  not  occur  even  at  pressures  of  1000  psi,  and  only  a 
small  amount  of  charring  was  observed  for  any  sample. 

It  was  concluded  that  the  use  of  ammonium  sulfate  as  a  total  replacement 
for  the  coarse-phase  perchlorate  in  RDS  501  yields  a  safe  propellant 
simulant . 


Duplication  of  Properties 

By  mutual  agreement  with  representatives  of  the  Air  Force,  the  solid  pro¬ 
pellant  parameters  to  be  duplicated  were: 

1.  Mechanical  Properties 

2.  Density 

3.  Thermal  Properties 

a.  specific  heat 

b.  coefficient  of  thermal  expansion 

c.  coefficient  of  thermal  conductivity 
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Hot  Wire  Ignition  Tests  of  Simulated  Propellant 
Containing  Only  Ammonium  Sulfate  in  the  Coarse- 
Oxidizer  Phase 


1252-5/11/62-SI 

Figure  77*  Gas  Generator  Ignition  Tests  on  Simulated  Propellants  Containing 
Only  Ammonium  Sulfate  in  the  Coarse  Oxidizer  Phase 


TABLE  18 


IGNITION  STUDIES  ON  RDS  501 :  THE  EFFECT  OF  REPLACEMENT 
OF  COARSE  PERCHLORATE  BY  AMMONIUM  SULFATE 


Percent  Sulfate  in 
Coarse  Phase 

Hot  ¥ire  Ignition 

Gas  Generator  Ignition 

Ambient 

500  psi 

1000  psi 

Ambient 

500  psi 

1000  psi 

0  (RDS  501) 

yes 

yes 

yes 

yes 

yes 

yes 

50 

no 

yes 

yes 

— 

— 

— 

75 

no 

yes 

yes 

— 

— 

— 

100 

no 

no 

no 

na 

no 

no 
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Mechanical  Properties.  The  mechanical  properties  of  live  RDS  501  and 
simulant  compositions  in  which  there  was  partial  or  complete  replacement 
of  the  coarse  perchlorate  by  sulfate  were  compared  by  use  of  the  JANAF 
tensile  properties  test  (Ref.  68) (  Table  19).  For  the  total  replacement 
simulant  and  RDS  501,  two  curing  schedules  were  employed  to  determine 
whether  time  or  temperature  affected  the  ability  to  duplicate  propellant 
properties.  A  very  satisfactory  match  of  properties  has  been  achieved. 
Table  20  compares  the  mechanical  behavior  over  a  temperature  range  of 
“70  to  4-170  F.for  the  live  propellant  and  the  simulant  containing  only 
ammonium  sulfate  in  the  coarse-oxidizer  phase.  Here  again,  the  duplica¬ 
tion  of  properties  is  demonstrated. 

Prior  work  at  Rocketdyne  had  shown  that  replacement  of  perchlorate  in 
both  the  coarse-  and  fine-oxidizer  phases  yields  a  product  which  does 
not  duplicate  the  mechanical  behavior.  Considering  the  results  described 
above,  it  is  felt  that  the  need  for  retention  of  perchlorate  in  the  fine- 
oxidizer  phase  has  been  established,  and  that  the  great  influence  of  this 
phase  on  the  propellant  mechanical  properties  demonstrated. 


Viscosity  of  Uncured  Propellant.  Although  propellant  viscosity  was 
not  designated  as  a  property  to  be  duplicated,  the  rheological  proper¬ 
ties  of  the  simulant  composition  become  important  if  it  is  to  be  used 
in  loading  motors  for  full-scale  testing.  The  viscosities  (at  150  F)  of 
RDS  501  and  several  simulants  are  compared  in  Table  21.  Castability 
(which  is  related  to  propellant  viscosity)  of  the  total  replacement  sim¬ 
ulant  can  be  made  to  match  RDS  501  if  the  former  is  cast  at  approximately 
190  F  (Table  22  ).  Since  there  is  no  safety  problem  with  the  simulant,  it 
appears  that  casting  at  this  temperature  is  feasible,  and  that  motors  can 
be  loaded  easily  with  simulant  propellant. 
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TABLE  19 


COMPARISON  OF  THE  MECHANICAL  PROPERTIES 
OF  RDS  501  AND  SIMULANTS 


Test 

0  percent  AS 
In  Coarse 
Phase 
(RDS  501) 

50  percent  AS 
In  Coarse 
Phase 

75  percent  AS 
In  Coarse 
Phase 

100  percent  AS 
In  Coarse 
Phase 

Young's 
Modulus,  psi 

1042* 

1225** 

1251* 

1432* 

1468* 

1009** 

Stress , 
maximum,  psi 

123.9* 

137.5** 

153* 

160* 

157.1* 

111.5** 

Elongation, 
maximum 
stress , 
percent 

20.8* 

20.3** 

19.0* 

20.0* 

21.1* 

22.3** 

Elongation 
at  break, 
percent 

21.3* 

21.2** 

10.5* 

20.5* 

21.6* 

23.4** 

*Curing  Schedule:  96  hours  at  170  F 

**Curing  Schedule:  72  hours  at  160  F 
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TABLE  20 


MECHANICAL  PROPERTIES  FOR  RDS  501  PROPELLANT 
AND  SIMULATED  RDS  501  PROPELLANT 


Property 

Temperature , 

F 

RDS  501 
Propellant 

RDS  501 
Simulant 

Maximum  stress  (S  ),  psi 

-70 

418. 6b 

521. 6C 

m 

+70 

137-5d 

111.5d 

+170 

112. 9C 

106. 5e 

Elongation  at  maximum  stress 

25. 2b 

(e  ),  percent 

-70 

24. 0C 

m 

+70 

22.3d 

+170 

13.7° 

14. 4e 

Elongation  to  break  (e^)j 

percent 

-70 

27. 4b 

25. 4C 

+70 

21. 2a 

23. 4d 

+170 

14. 3C 

15.  le 

Modulus  (initial  tangent), 

7,l63b 

i 

psi 

-70 

10,003c 

+70 

l,225d 

1 , 009d 

+170 

1 ,039c 

1 , 029e 

Notes : 

aTotal  replacement  of  perchlorate  by  ammonium  sulfate  in  the  coarse 
oxidizer  phase 


Average 

of 

22 

CAverage 

of 

20 

^Average 

of 

10 

e Average 

of 

19 
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TABLE  21 


VISCOSITY  OF  RDS  501:  EFFECT  OF  AMMONIUM 
SULFATE  IN  COARSE  OXIDIZER  PHASE 


Percent  Sulfate  in  Coarse  Phase 

Viscosity  at  150  F,  cps 

0 

975,000 

50 

1,067,000 

75 

1,125,000 

100 

1,250,000 

TABLE  22 

EFFECT  OF  TEMPERATURE  ON  VISCOSITY  OF 
RDS  501  SIMULANT* 


Temperature,  F 

Viscosity,  cps 

150 

1,250,000 

187 

800,000 

210 

500,000 

*Only  ammonium  sulfate  present  in  the  coarse  phase. 


147 


Density.  Density  measurements  were  performed  using  the  basic  principle 
of  Archimedes  with  2,  2,  4-trimethylpentane  as  the  immersion  fluid.  Pro¬ 
pellant  ingredients  are  insoluble  in  this  fluid.  All  density  tests  were 
performed  at  77  F.  The  density  of  RDS  501  was  determined  to  be  1.75  g/cc; 
the  density  of  the  simulant  was  1.67  g/cc. 


Replacement  of  Coarse  Perchlorate  by  Other  Ingredients.  The  density 
of  the  simulant  containing  ammonium  sulfate  in  the  coarse  phase  is  5  per¬ 
cent  lower  than  RDS  501.  For  this  reason,  an  investigation  on  the  use  of 
more  dense  ingredients  was  begun.  Potassium  chloride  (l.98  g/cc),  sodium 
chloride  (2.17  g/cc),  and  ammonium  chloride  ( 1 . 94  g/cc)  were  selected. 
These  salts  were  ground,  screened,  and  dried  so  that  no  aberration  in 
results  could  be  attributed  to  moisture  or  the  introduction  of  these 
ingredients  into  the  fine  phase. 

Propellant  containing  ammonium  chloride  could  not  be  made  at  the  desired 
solids  loading  level.  Castability  difficulties  were  encountered  even  at 
much  lower  solids  loading.  The  potassium  chloride  propellant  did  not 
cure  at  all.  Only  the  propellant  containing  sodium  chloride  gave  a  satis¬ 
factory  cure.  The  mechanical  properties  (at  70  F)  of  this  simulant,  the 
ammonium  sulfate  version,  and  RDS  501  are  compared  in  Table  23~  The 
sodium  chloride  simulant  does  not  duplicate  RDS  501  as  well  as  the  ammon¬ 
ium  sulfate  composition.  For  this  reason,  this  approach  was  not  pursued 
in  evaluating  how  well  it  matches  the  thermal  properties  of  RDS  501. 

To  determine  the  curing  reaction  of  binders  in  the  presence  of  different 
fillers,  a  study  was  made  of  the  absorption  of  the  MAP0  crosslinker  by 
ammonium  sulfate,  potassium  chloride  and  sodium  chloride.  The  latter 
two  salts  absorb  the  crosslinker  to  an  extent  three  times  greater  than 
the  sulfate.  Therefore,  it  was  not  surprising  that  the  propellant 
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TABLE  23 


COMPARISON  OF  MECHANICAL  PROPERTIES 
OF  RDS  501  AND  SIMULANTS 


Property 

Temperature , 

F 

RDS  501 
Propellant 

RDS  501 
Simulant  * 

RDS  501 

2 

Simulant 

Maximum  Stress 
(Sm),  pai 

+  70 

137- 53 

m.53 

104.0 

Elongation  at 

Maximum  Stress 
(e^),  percent 

+  70 

20. 33 

22. 33 

14.1 

Elongation  to 

Break  (e^) ,  percent 

+  70 

21. 23 

23. 43 

18.8 

Modulus 

(initial  tangent),  psi 

+  70 

3 

1,225 

1,0093 

1,069 

Notes : 

1.  Total  replacement  of  perchlorate  in  the  coarse  oxidizer  phase  by 
ammonium  sulfate 

2.  Total  replacement  of  perchlorate  in  the  coarse  oxidizer  phase  by 
sodium  chloride 

3.  Average  of  10  values 
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containing  the  potassium  chloride  did  not  cure.  The  curing  of  the  sodium 
chloride  simulant,  which  may  be  only  a  partial  cure  (note  the  poorer 
mechanical  properties)  may  be  due  to  the  ability  of  the  polymer  to  desorb 
the  crosslinker  more  readily  from  sodium  chloride  than  from  potassium 
chloride.  The  experimental  difficulties  involved  in  proving  this  point 
precluded  any  attempt  under  this  contract.  This  point,  while  remaining 
a  conjecture,  appears  to  be  a  reasonable  explanation  of  the  results. 


Thermal  Properties.  The  thermal  properties  of  solid  propellants  are  not 
measured  routinely.  Thus,  this  phase  of  the  program  not  only  involved 
the  measurement  of  such  parameters,  but  also  required  an  evaluation  of 
the  applicability  of  the  methods  which  had  been  developed  for  nonpropel¬ 
lant  ingredients. 


Thermal  Expansion.  The  standard  ASTM  method  of  determining  the 
coefficient  of  linear  thermal  expansion  was  used  for  evaluation.  This 
method  employs  a  vycor-tube  dilatometer  apparatus.  The  coefficients  for 
copper  and  plexiglass  were  determined  to  verify  the  test  technique. 

Initial  measurements  at  elevated  temperatures  were  made  using  a  tempera¬ 
ture-controlled  bath  with  later  tests  conducted  in  an  environmental 
chamber.  This  latter  apparatus  provided  greater  accuracy  by  permitting 
continuous  measurements  from  -70  F  to  +  170  F,  thus  reducing  the  effect 
of  the  small  static  errors  inherent  in  the  system.  Propellant  samples 
of  approximately  0.4  inch  in  diameter  by  4.5  inches  long  were  conditioned 
at  each  temperature  step  for  at  least  1  hour.  The  resulting  displacement 
was  monitored  until  constant  dilatometer  readings  were  obtained.  Each 
test  consisted  of  a  minimum  of  1  cycle  through  the  entire  temperature 
range,  with  additional  measurements  conducted  in  the  vicinity  of  transi¬ 
tion  areas. 
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From  these  measurements,  the  expansion  coefficients  «  were  calculated 
via  Eq.  49: 


«  -  urmj  (M> 

The  results  are  described  in  Table  24  and  Fig.  78.  The  simulant  coef¬ 
ficient  differs  from  that  of  the  live  propellant  by  approximately  20  per¬ 
cent.  In  addition,  while  the  coefficient  for  the  live  EDS  501  was  linear 
over  the  entire  temperature  range  (from  -70  to  +170  F)  the  simulant  expan¬ 
sion  coefficient  exhibited  a  transition  at  -32  F.  Thus,  for  the  simulant, 
Of  decreases  from  3.62  x  10  ^  (-32  to  +170  F)  to  1.26  x  10  ^  (-32  to  -70  F). 
A  similar  transition  point  occurred  with  a  simulated  polyurethane  propel¬ 
lant.  Although  both  polyurethane  formulations  exhibited  higher  coeffici¬ 
ents  over  the  major  portion  of  the  temperature  range,  the  simulated  poly- 

-5 

urethane  propellant  coefficient  changed  radically  from  6.55  x  10  to 
0.9  X  10“5  at  -42  F. 


Specific  Heat.  The  test  method  for  the  determination  of  specific 
heat  was  a  slightly  modified  version  of  the  ASTM  procedure,  which  is 
based  on  the  classical  method  of  mixtures.  The  modification  employed 
was  the  use  of  a  threaded  aluminum  capsule  (Fig.  79)  with  a  thermocouple 
well,  rather  than  the  standard  thin-wall  brass  container.  This  change 
was  made  for  economy  because  the  propellant  was  cast  into  the  capsule, 
and  the  capsule  was  considered  expendable  after  one  test.  In  the  test, 
a  temperature-conditioned  specimen  of  known  mass  is  immersed  in  a  water- 
filled  calorimeter  at  a  lower  temperature  than  that  of  the  propellant. 

The  resulting  temperature  rise  to  thermal  equilibrium  and  the  specimen 
weight  are  used  to  compute  specific  heat.  Evaluation  of  the  test  appa¬ 
ratus  and  technique  was  conducted  over  a  temperature  range  of  75  to  170  F 
using  both  copper  and  plexiglass. 
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TABLE  24 


COEFFICIENT  OF  LINEAR  EXPANSION 
-70  to  170  F 


Materia  1 

/ ttl 

in/in-F 

Temperature 
Range,  F 

Transition 
Temperature, F 

*2 

in/ in-F 

RDS  501  Propellant  & 

4.  55xio”5 

-70  to  170 

None  >  -70 

— 

RDS  501  Simulant  b  C 

3.62x10-5 

-32  to  170 

-32 

1.26xl0~5 

Polyurethane  Propellant  9 

7.59x10-5 

-48  to  170 

-48 

4.62xl0~5 

Polyurethane  Simulant 

6.55x10-5 

-42  to  170 

-42 

0.9x10" 5 

Notes : 

a.  Average  of  five  tests  from  tvo  batches,  cr  ±  5  percent 

b.  Total  replacement  of  perchlorate  by  ammonium  sulfate  in  the  coarse  phase 

c.  Average  of  four  tests  from  three  batches,  O'  ±  7  percent 

d.  Single  test 
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70 


170 


"*u  TEMPERATURE,  F 

♦Total  Replacement  of  Coarse-Phase  Perchlorate  by 
Ammon ium  Sulfate 

Figure  78.  Typical  Propellant  Thermal  Expansions, 
-70  to  170  F 


153 


The  average  specific  heat  of  several  propellants  and  propellant  simulants 
was  determined  for  the  temperature  range  of  66  to  170  F.  Table  25  sum¬ 
marizes  these  results,  which  show  that  the  oxidizer  (the  major  portion  of 
the  propellant)  is  the  major  factor  in  determining  the  specific  heat. 

Thus,  replacement  of  ammonium  perchlorate  by  ammonium  sulfate  in  the 
coarse  phase  results  in  an  increase  in  specific  heat.  The  specific  heat  of 
the  simulant,  therefore,  deviates  by  approximately  20  percent  from  that  of 
the  live  propellant. 

To  check  these  results  and  to  determine  the  specific  heat  of  propellant 
without  the  hazards  involved  in  making  measurements  of  a  live  propellant 
the  specific  heat  was  calculated  from  the  specific  heats  of  the  propel¬ 
lant  ingredients.  As  it  can  be  seen  from  Table  25  excellent  agreement 
between  the  experimental  and  calculated  values  was  obtained.  This  non- 
hazardous  method  for  determining  propellant  specific  heat  appears  feasible. 


Thermal  Conductivity.  The  use  of  standard  methods  for  determining 
thermal  conductivity  or  thermal  diffusivity  has  not  been  widespread; 
consequently  many  of  the  methods  used  or  proposed  are  not  particularly 
applicable  to  solid  propellants.  The  principal  objective  of  this  phase 
of  the  program,  therefore,  was  to  evaluate  existing  methods  critically 
for  accuracy  and  reliability,  and,  if  possible,  to  develop  a  technique 
which  would  circumvent  the  hazards  and  inadequacies  inherent  to  most  of 
these . 

Since  steady-state  as  well  as  transient  heat  properties  are  important, 
the  major  effort  concerned  the  measurement  of  thermal  diffusivity. 

This  nonsteady-state  property,  together  with  specific  heat  and  density 
values,  permits  the  calculation  of  thermal  conductivity. 
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TABLE  25 


A'VERAGE  SPECIFIC  HEAT,  -66  TO  170  F, 
OF  EDS  501  AND  SIMULANT 


Material 

No.  of 
Batches 

No.  of 
Tests 

^P7 

Btu/lb-F 

EDS  501  Propellant 

3 

7 

0.258  ±10  percent 

EDS  501  Simulant* 

6 

15 

0.309  ±10  percent 

EDS  501  Binder 

1 

3 

0.274  ±10  percent 

Polyurethane  Propellant 

1 

2 

0.285  ±2 

Polyurethane  Simulant* 

1 

2 

0.325  ±5 

Ammonium  Perchlorate) 

0.268 

Ammonium  Sulfate  )  Literature  Values 

0.343 

Aluminum  ) 

0.218 

EDS  501  Propellant)  Calcula 

1 

ted  from  ingredient 

0.261 

EDS  501  Simulant*  )  mass  fractions  and  literature 

0.298  1 

values 

* 

Total  replacement  of  perchlorate  in  the  coarse-oxidi zer  phase 
by  ammonium  sulfate. 
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The  specific  techniques  investigated  during  this  program  were  the  methods 
of  Chung  and  Jackson  (Ref.  69),  Schneider  (Ref.  70  ),  and  Brown  and 
Marco  (Ref.  71  )•  Several  additional  methods  were  reviewed  briefly,  but 
were  rejected  as  not  being  suitable  for  propellants. 


One  or  more  of  the  following  factors  is  inherent  to  the  unsatisfactory 
methods  and  is  sufficient  justification  for  not  considering  these  tech¬ 
niques.  These  factors  are: 


1.  Exotic  specimen  geometry;  i.e.,  hollow  sphere 

2.  Costly  hardware  or  instrumentation 

3.  Gross  temperature  errors 

4.  Invalid  assumptions 

5.  Test  conditions  not  suitable  for  propellants 


Basically,  the  Chung  and  Jackson  method  differs  from  the  other  two  in 
the  geometry  of  the  test  specimen  (cylindrical  vs  cubic).  In  this 
method,  the  specimen  is  brought  to  a  constant  initial  temperature  in  a 
vacuum  oven.  Upon  reaching  this  predetermined  temperature,  the  specimen 
is  immediately  transferred  to  a  heat  transfer  cell  where  water  is  allowed 
to  flow  by  at  a  prescribed  flowrate.  The  specimen  should  have  a  length- 
to-diameter  ratio  of  8  or  more.  A  thermocouple  is  placed  anywhere  along 
the  centerline  of  the  cylinder  because  it  is  assumed  that  only  radial 
heat  transfer  occurs.  This  is  a  safe  assumption  only  when  the  length-to- 
diameter  ratio  is  8  or  more. 


The  other  methods  are  similar  in  procedure  in  that  the  sample  is  brought 
to  an  initial  constant  temperature,  and  then  immersed  in  a  liquid  quenching 
bath.  The  methods  differ  primarily  in  the  treatment  of  the  data. 
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The  initial  evaluation  of  the  three  methods  was  conducted  using  inert 
materials  having  thermal  properties  very  similar  to  polybutadiene-type 
propellants.  Three  materials,  Teflon,  Kel-F,  and  plexiglass  were  tested 
in  both  cylindrical  and  cubical  configurations.  Adequate  sealing  of  the 
first  two  materials  to  a  centerline  thermocouple  was  not  possible.  Small 
insulating  air  pockets  in  the  drilled  thermocouple  wells  could  not  be 
satisfactorily  eliminated.  These  two  materials  were  abandoned  as  stand¬ 
ards  because  of  the  uncertainty  of  temperature  measurements  under  these 
conditions.  The  sealing  of  the  thermocouples  in  plexiglass  was  accom¬ 
plished  using  a  slush  of  plexiglass  shavings  and  methylene  chloride,  and 
resulted  in  a  firm  homogeneous  seal. 

The  standardization  of  the  Chung  and  Jackson  apparatus  (Fig.  80 )  was 
carried  out  using  a  precisely  machined  plexiglass  cylinder  sealed  to 
the  cover  plate  of  the  heat  transfer  cell  with  an  epoxy  resin.  Length- 
to-diameter  ratios  ranging  from  6  to  8  were  tested  at  temperatures  up 
to  170  F. 

Both  the  method  of  Schneider  and  that  of  Brown  and  Marco  require  that 
a  sample  be  brought  to  constant  initial  temperature  and  subsequently 
quenched  in  a  cooler,  constant-temperature  bath.  It  is  evident  that  a 
large  temperature  difference  is  necessary  to  ensure  accurate  results. 

Due  to  the  relatively  low  thermal  diffusivity  of  propellants,  several 
minutes  of  exposure  to  the  cooler  environment  were  required  before  the 
heat  flux  boundary  reached  the  thermocouple.  Also,  it  was  observed  that 
for  accurate  results  it  was  not  necessary  that  the  specimen  be  maintained 
in  the  quenching  bath  until  thermal  equilibrium  was  established.  There¬ 
fore,  in  the  standard  test,  3  minutes  were  allowed  for  the  heat  flux 
boundary  to  match  the  thermocouple  before  measurements  were  begun,  and 
after  13  more  minutes  no  further  temperature  recordings  were  taken. 
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Figure  80.  Chung  and  Jackson  Apparatus 


The  proper  specimen  dimensions  were  determined  by  successive  testing  of 
0.7,  1.0,  and  2.0-inch  plexiglass  cubes.  Testing  shoved  that  a  2.0-inch 
cube  was  necessary  to  ensure  accuracy  vithin  97*5  percent  (Fig.  81 ) .  A 
1.0-inch  cube  was  only  75  to  80  percent  accurate,  vhile  the  0.7-inch 
cube  was  100  percent  in  error.  The  smaller  samples  were  tested  first  to 
decrease  operating  hazards  by  minimizing  the  amount  of  propellant  present. 

Because  specimen  preparation  is  a  major  source  of  error,  a  detailed  pro¬ 
cedure  vas  evolved  as  follows: 

1.  Ingredients  are  mixed  at  elevated  temperatures,  cast  into 
block  configurations;  then  cured  for  96  hours  at  170  F. 

These  batches  were  cast  with  30-gage  iron-constantan  thermo¬ 
couples  in  place. 

2.  All  cured  castings  were  X-rayed.  A  minimum  of  two  views  was 
obtained  for  each  casting. 

3.  X-ray  plates  were  examined  for  void  indications  and  for  accur¬ 
acy  of  thermocouple  placement.  Castings  with  obvious  faults 
were  discarded. 

4.  Each  satisfactory  casting  was  then  carefully  machined  to  the 
dimensions  required. 

5.  Test  specimens  were  stored  at  77  F  until  used. 


Comparison  of  the  Three  Test  Methods.  In  the  application  of  the 
Chung  and  Jackson  method,  several  conditions  were  encountered  which, 
regardless  of  the  specimen  length-t o-diaraeter  ratio,  produced  cumula¬ 
tive  errors  in  the  measurements  as  high  as  40  percent.  For  example, 
in  this  technique,  the  conditioned  specimen  must  be  placed  into  the  heat 
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PLEXIGLASS  CUBE  EDGE,  INCHES 


3 


2 


1 


20  40  60  80  100 

PERCENT  DEVIATION  OF  THERMAL  DIFPUSIVITY 
FROM  LITERATURE  VALUE 


Figure  81.  Variation  of  Plexiglass  Thermal  Diffusivity 
Measurement  vith  Specimen  Size 
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transfer  cell  and  securely  fastened  into  place.  This  presented  a  prob¬ 
lem  in  the  application  of  a  correction  factor  to  compensate  for  the  time 
lag  and  heat  loss  occurring  before  the  recording  of  the  cooling  curve 
could  be  started.  The  other  two  methods  use  quenching,  and  do  not  suf¬ 
fer  from  this  handicap. 

In  general,  cylindrical  samples  present  more  difficulty  in  the  correla¬ 
tion  of  experimental  data.  This  is  attributed  to  the  fact  that  various 
methods  of  correlation  depend  upon  the  assumption  of  an  infinite  or 
semi-infinite  cylinder.  The  distinction  between  the  two  is  not  clearly 
defined.  One  other  obvious  disadvantage  inherent  to  cylindrical  speci¬ 
mens  having  the  high  length-to-diameter  ratios  mandatory  for  this  method 
is  the  additional  mass  of  propellant  required  to  achieve  any  increase  in 
the  length  of  the  radial  heat. flow  path.  This  method  was  abandoned 
after  preliminary  standardization  tests. 

A  study  of  all  the  data  obtained  indicated  that  Schneider1 s,  and  Brown 
and  Marco’s  methods  suffer  from  a  time-dependency  factor.  This  factor 
is  not  too  large  in  Schneider’s  method  but  is  more  pronounced  in  Brown 
and  Marco’s  method.  Several  minutes  must  be  allowed  in  both  methods 
before  an  accurate  calculation  can  be  made.  For  a  temperature  range  of 
150  F,  a  6-  to  8-minute  period  following  the  initial  3~  to  4-minute  delay 
will  yield  good  results.  After  10  minutes  total  elapsed  time,  Brown  and 
Marco’s  correlation  charts  become  asymptotic  and  cannot  be  read  accur¬ 
ately.  Schneider’s  method  does  not  deviate  as  rapidly  and  gives  accurate 
results  for  approximately  5  additional  minutes.  A  comparison  is  shown 
in  Fig.  82. 

The  cubic  samples  used  for  these  tests  are  a  definite  advantage.  They 
can  be  analyzed  by  methods  applicable  to  any  rectangular  parallelopiped . 
In  addition,  a  cube  allows  several  simplifying  assumptions.  Larger  size 
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Figure  82.  Comparison  of  Brown  &  Marco  Method  vs  Schneider  Method  for  1-inch  Plexiglass  Tube 


specimens  offer  a  definite  advantage  in  increased  accuracy.  Jakob 
(Ref.  72)  has  pointed  out  the  errors  introduced  in  traversing  a  range 
of  dimensions. 

The  Schneider  method,  therefore,  clearly  demonstrated  superior  accuracy 
and  reproducibility  over  the  other  two  candidates.  All  propellant  ther¬ 
mal  diffusivity  tests  were  conducted  using  Schneider* s  method.  Possible 
sources  of  error  in  the  method  are:  (l)  voids,  (2)  leaching  of  soluble 
salts  by  water,  and  (3)  invalid  temperature  measurements. 

The  undesirable  effects  of  voids  can  never  be  completely  eliminated, 
but  this  can  be  held  to  a  minimum  by  careful  processing  of  test  speci¬ 
mens.  Some  authors  have  contended  that  a  0.1-percent  difference  in  void 
content  would  change  the  thermal  diffusivity  by  10  percent.  This  appears 
to  be  too  high  an  estimate. 

Nevertheless,  voids  would  act  as  a  thermal  barrier  and  their  effect  upon 
the  thermal  diffusivity  must  be  taken  into  account.  By  performing  bulk 
modulus  tests  on  a  sample  from  each  batch  of  propellant,  one  could  gain 
some  insight  into  the  true  variation  of  thermal  diffusivity  with  void 
content . 

Leaching  of  ammonium  perchlorate  was  observed  upon  immersion  of  the  pro¬ 
pellant  samples  in  the  water  baths.  Investigation  showed  that  this  was 
a  surface  effect  only  and  that  0.25  grams  of  ammonium  perchlorate  was 
leached  out  of  a  220-gram  sample.  This  amount  is  not  large  enough  to 
alter  the  test  results.  Cold  tests  using  Freon  and  liquid  nitrogen  did 
not  present  any  obvious  problems. 


The  use  of  small  thermocouples  (30-gage  iron-constantan)  cast  in  posi¬ 
tion,  together  with  a  precision  potentiometer  for  measuring  emf ,  mini¬ 
mize  the  possibility  of  invalid  temperature  measurements.  Smaller 
thermocouple  wire  with  the  lead-in  wire  coiled  in  the  propellant  would 
further  increase  the  accuracy  of  temperature  measurements. 

The  thermal  dif fusivities  of  HDS  501  and  the  simulant  in  which  all  the 
coarse  perchlorate  was  replaced  by  sulfate  were  0.00971  and  0.00856  sq 
ft/hr,  respectively,  over  the  temperature  range  of  77  to  200  F.  Using 
these  values,  and  the  densities  and  specific  heats  reported  above,  the 
thermal  conductivities  were  calculated  (0.273  and  0.274  Btu/ft  F). 

The  thermal  diffusivity  of  the  simulant  was  measured  down  to  -100  F, 
and  is  very  nearly  linear  throughout  (Fig.  83).  The  equation  for  the 
simulant  diffusivity  is: 

«  =  1.067  x  10-2  -  1.325  x  10-5  T 

where 

T  =  degrees,  F 

O'  =  sq  ft/hr 

This  equation  should  serve  as  an  excellent  approximation  for  the  live 
propellant  as  well. 


Table  26  is  a  summary  of  all  the  thermal  properties  data  accumulated  for 
the  live  and  simulated  propellant.  It  is  felt  that  a  good  duplication 
of  propellant  parameters  has  been  obtained,  particularly  in  relation  to 
thermal  diffusivity.  Thermal  diffusivity  is  probably  of  major  importance 
in  escape  or  reentry  situations  where  transient  heat  transfer  conditions 
prevail . 
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Simulated  Propel 
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